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Abstract. We review the theory of the temperature anisotropy and polarization of the cosmic 
microwave background (CMB) radiation, and describe what we have learned from current CMB 
observations. In particular, we discuss how the CMB is being used to provide precise measurements 
of the composition and geometry of the observable universe, and to constrain the physics of the early 
universe. We also briefly review the physics of the small-scale CMB fluctuations generated during 
and after the epoch of reionization, and which are the target of a new breed of arcminute-resolution 
instruments. 
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1. INTRODUCTION 

The last two decades have seen remarkable advances in observational cosmology and 
quantitative cosmological constraints now come from a number of complementary 
probes. Amongst these, the cosmic microwave background temperature and polarization 
anisotropics have played a particularly important role. The physics is very well under- 
stood - the primary contribution to the anisotropies is accurately described by linear 
perturbation theory - and so theoretical predictions are believed to be accurate to better 
than 1%. Furthermore, the CMB currently provides the most direct observational link to 
the physics of the early universe since the fluctuations were mostly imprinted around the 
time of recombination. 

The physics of the CMB anisotropies is now a textbook subject [1]. The basic picture 
is that the angular variations in temperature that we observe today are a snapshot of the 
local properties (density, peculiar velocity and the total gravitational potential) of the 
gas of CMB photons at the time the primordial plasma recombined at redshift z ^ 1 100. 
These quantities are related to the primordial perturbations, plausibly generated during 
an epoch of cosmological inflation, by the rather simple acoustic physics of the pre- 
recombination plasma. The CMB anisotropies therefore encode information on the 
primordial perturbation itself, as well as the matter composition and geometry of the 
universe. By mapping the microwave sky, cosmologists thus hope to answer some of the 
biggest questions in physics. 

Large-angle temperature anisotropies were first discovered by COBE-DMR in 
1992 [2]. Since then, tens of instruments have mapped the anisotropies and we now 
have high quality measurements of the statistics of the anisotropies over three decades 
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FIGURE 1. Fluctuations in the CMB temperature, as determined from five years of WMAP data, about 
the average temperature of 2.725 K. Reproduced from [3]. 



of angular scales. The current state-of-the-art dataset is provided by COBE's successor, 
the WMAP satellite, which released its five-year results in March 2008 [3]. A map of the 
temperature fluctuations from WMAP, after removing the effect of microwave emission 
from our galaxy, is shown in Fig. 1 The WMAP results are beautifully fit by a simple 
cosmological model with flat space, baryons and cold dark matter making up 4% and 
21% respectively of the energy density with the remaining 75% in some form of dark 
energy, and Gaussian-distributed, adiabatic primordial density perturbations with an 
almost, but not quite, scale-invariant spectrum. The accuracy of these numbers is quite 
remarkable by cosmological standards; for example, the baryon density is measured 
to better than 3% and space is flat to better than 1%. Inflation continues to stand up 
to exacting comparison with observations. Indeed, the quality of the data is now good 
enough to rule out several popular inflation models from the plethora that theorists have 
proposed. The community now eagerly awaits the launch of Planck, the third-generation 
CMB satellite, in mid-2009. This will image the sky in nine wavelength bands with 
a best resolution below 5 arcmin. Forecasts indicate that Planck should measure sev- 
eral important parameters to better than 1%, and place tight constraints on inflation 
models [4]. 

The CMB temperature fluctuations are complemented by its linear polarization [5] . 
Polarization is generated by scattering of anisotropic radiation and for the CMB this oc- 
curs at two epochs: around recombination and, later, from the epoch of reionization [6]. 
The small polarization signal (r.m.s. ~ 5 jlK) was first detected in 2002 by the DASI in- 
strument [7], and has subsequently been measured by several groups. The signal-to-noise 
of these measurements is still poor compared to the temperature anisotropics, but already 
important constraints have emerged on the epoch of reionization from WMAP [8, 9]. 
Polarization also provides the only viable route in the near-future for detecting the back- 
ground of gravitational waves predicted from inflation in models with large, i.e. super- 
Planckian, field values [10, 11]. A positive detection would be particularly interesting 
for attempts to realise inflation in string theory, and would rule out some alternatives 
to inflation such as the cyclic model [12]. Since Planck may not have the sensitivity to 
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make a definitive test of gravitational waves, a new generation of ambitious sub-orbital 
instruments that will deploy hundreds or even thousands of superconducting detectors 
to beat down the noise are under construction. These exciting instruments aim to probe 
the highly uncertain physics of energies near 10^^ GeV. 

In addition to these primary anisotropics, there are a number of secondary processes 
that add further structure to the CMB at late times (see Ref. [13] for a recent review). 
These include gravitational redshifts induced by the linear evolution of the gravitational 
potentials once dark energy becomes dynamically important [14], gravitational lensing 
of CMB photons by the large-scale distribution of matter [15, 16, 17], and various scat- 
tering effects that become operative after the onset of reionization [18, 19, 20]. While the 
observation and interpretation of the primary anisotropics is now a mature field, for the 
secondary processes, in which the CMB is used as a back-light, this is just beginning. A 
new breed of instruments with arcminute resolution are now surveying these secondary 
effects with the hope of improving our understanding of the epoch of reionization and 
the growth of structure over a wide range of redshifts. If these new observations are to 
fulfill their scientific potential, accurate modelling of secondary anisotropics is required. 
This properly requires large numerical simulations and is currently a very active area of 
CMB theory. 

In these lecture notes we outline the basic physics of the primary temperature and 
polarization anisotropics of the CMB, and discuss how observations of these are being 
used to constrain cosmological models and the physics of the early universe. For earlier 
reviews that are similar in spirit, see for example [21, 22, 23]. We also briefly review 
secondary anisotropics and look forward to what we may learn from the new high- 
resolution surveys. 

2. INTRODUCTION TO CIMB TEIMPERATURE ANISOTROPIES 

The CMB radiation has an exquisite blackbody (thermal) spectrum, with very nearly 
the same temperature of 2.725 K in all directions on the sky [24, 25]. The largest 
temperature variation is a dipole due to our motion relative to the CMB. Initially, ignore 
the small effect of the temperature fluctuations due to the primordial perturbations. Then 
the CMB is isotropic in its rest frame with a Lorentz-invariant distribution function^ 
/(pM) oc i/(e£cMB/7cMB _ 1) where Eqmr is the photon energy in the CMB frame and 
TcMB = r is the isotropic temperature. If we observe a photon with energy E and 
direction e, then the relativistic Doppler shift gives £'cmb = Ey{ \ -|- v • e), where v is 
our velocity relative to the CMB. In terms of the energy and direction in our reference 
frame, 

f^P"") " ,£y(i+v.e)/rcMB _ 1 • 

This still looks like a blackbody along any direction but the observed temperature varies 
over the sky as r(e) — 7cmb/[7(1 + v • e)] pa 2cmb(1 — v • e) for |v| -C 1, i.e. a dipole 
anisotropy. The measured dipole [26] implies the solar system barycentre has speed 



The one-particle distribution function gives the number of photons per proper phase-space volume. 
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3.68 X 10^ ms^ relative to the CMB. It is clear from Eq. (1) that relative motion also 
produces quadrupole anisotropies at 0(|vp) and CMB data is routinely corrected for 

this (small) effect. Reinstating the contributions of the primordial perturbation to the 
temperature anisotropies, these are at the level of one part in 10^ and dominate over the 
relative- velocity effects for the observed quadrupole and higher multipoles. 



2.1. CMB observables 

It is convenient to expand the temperature fluctuation in spherical harmonics, 

Ar(n)/f = J^M«(n), (2) 

Im 

with — ( — l)'"r/_„ since the temperature is a real field. The sum in Eq. (2) runs 
over / > 1, but the dipole (/ = 1) is usually removed explicitly when analysing data 
since it depends linearly on the velocity of the observer. Multipoles at / encode angular 
information with characteristic scale (i.e. period) ~ 2k /I. 

The statistical properties of the fluctuations in a perturbed cosmology can be expected 
to respect the symmetries of the background model. In the case of Robertson-Walker 
models, this means the statistics should be independent of spatial position (homoge- 
neous) and invariant under rotations of the field. Under a rotation specified by Euler 
angles a, j8 and y, the T/^ transform to Lm'-DLm'("'J^' ^)^/'n' ^^^^^ ^ Wigner 

rotation matrix. Demanding invariance under rotations fixes the second-order statistics 
to be of the form 

{TlmTijm') = CfSii'dmm! , (3) 

which defines the power spectrum Cf . The angle brackets in this equation denote the 
average over an ensemble of realisations of the fluctuations. The simplest models of 
inflation predict that the fluctuations should also be Gaussian at early times, and this is 
preserved by linear evolution of the small fluctuations. Equation (3) tells us that the Tifn 
are uncorrelated for different / and m and, if Gaussian, they are also independent. In this 
case, the power spectrum provides the complete statistical description of the temperature 
anisotropies. For this reason, measuring the anisotropy power spectrum has, so far, been 
the main goal of observational CMB research. 

If we were given ideal (i.e. noise-free) measurements of the temperature over the full 
sky, the CMB power spectrum could be estimated by 

Cf-^B^lml'- (4) 

This is an unbiased estimator of the true ensemble Cf, but has an irremovable cosmic 
variance due to the finite number (2/ -|- 1) of modes that we can observe. If the temper- 
ature anisotropies are Gaussian distributed, Cf has a distribution with 2/ -|- 1 degrees 
of freedom and the cosmic variance is [27] 

var(Cn = ^(Cn^ (5) 



Lecture notes on the physics of cosmic microwave background anisotropies 



March 30, 2009 



4 



6000 









5000 
^ 4000 
3 3000 






■ WMAP 5-year 
• ACBAR 

■ BOOMERANG03 






t 


^~ 2000 






Mi 


1000 





















1 


10 100 


500 1000 1500 2000 2500 



FIGURE 2. The temperature anisotropy power spectrum, &i = l{l + \)Cj /{2n), as measured by 
WMAP5 (blue; [29]), ACBAR (black; [30]) and BOOMERanG (cyan; [31]). The line is the best-fit flat, 
ACDM model to the WMAP5 and ACBAR data. Reproduced from [30]. 

In practice, estimating the power spectrum is complicated by a number of real-world 
complexities such as partial sky coverage and instrumental noise. The critical review by 
Efstathiou [28] provides a useful entry point to the large literature on practical methods 
for estimating the temperature power spectrum. 

We show in Fig. 2 a selection of recent measurements of Cj taken from [30]. WIMAP 
and ACBAR are currently the most constraining of all measurements of Cj on large 
and small scales, respectively, but will be improved upon significantly by Planck for 
/ > 500. (W1VIAP5 is cosmic-variance limited to I = 530 [29].) Also shown in Fig. 2 
is the theoretical power spectrum for the best-fit flat, ACDIVI model to the W1VIAP5 
and ACBAR data. The agreement is striking, with five acoustic peaks clearly delimited. 
iVloreover, the agreement between different experiments (and technologies) is excellent 
where they overlap in / coverage. Understanding the main features in the Cf spectrum 
is one of the principal aims of these lecture notes. 

The correlation between the temperature anisotropics along two directions evaluates 

to 

(Ar(ni)Ar(n2)) = T^Y^^^Cf Pi {cos 0) , (6) 

which depends only on the angular separation 6 as required by rotational invariance. 
Here, Pi{x) are the Legendre polynomials. The mean- square temperature anisotropy is 

i^T') = T'l^^-^Cf ^ T'J din/ , (7) 

so that the quantity /(/ -|- l)Cf /2k, which is conventionally plotted, is approximately the 
power per decade in I of the temperature anisotropies. 

Higher-Order Statistics. In Gaussian theories, all higher-order correlations of the 
Ti,„ can be reduced to the power spectrum by Wick's theorem. In non-Gaussian theories. 
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the form of the correlators is restricted by rotational invariance. Here, we shall just 
consider the three-point function which must take the form 



where the right-hand side involves the Wigner 3 j symbol (which enforces mi+m2 + 
m3 = 0) and the angular bispectrum Bi^i^i^. If we further assume that the generation 
and propagation of cosmological perturbations respects parity symmetry, the bispectrum 
must then vanish for /i -|- /2 + h odd. Parity invariance further ensures that the bispectrum 
is real and invariant under arbitrary perturbations of its indices. There is, as yet, no strong 
evidence for a non-zero CMB bispectrum; see Sec. 5.3.3. For a discussion of the CMB 
four-point function, or trispectrum, see [32]. 



The high temperature of the early universe maintained a low equilibrium fraction of 
neutral atoms, and a correspondingly high number density of free electrons. Coulomb 
scattering between the ions and electrons kept them in local kinetic equilibrium, and 
Thomson scattering of photons tended to maintain the isotropy of the CMB in the baryon 
rest frame. As the universe expanded and cooled, the dominant element hydrogen started 
to recombine when the temperature fell below ~ 4000 K. This is a factor of 40 lower than 
might be anticipated from the 13.6-eV ionization potential of hydrogen, and is due to the 
large ratio of the number of photons to baryons. Through recombination, the mean-free 
path for Thomson scattering grew to the horizon size and CMB photons "decoupled" 
from matter. More precisely, the probability density that photons last scattered at some 
time defines the visibility function. This is now known to peak 380 kyr after the big 
bang [3], with a width ~ 120 kyr. Since then, CMB photons have propagated relatively 
unimpeded for 13.7 Gyr, covering a comoving distance %* = 14.1 Gpc. The distribution 
of their energies carries the imprint of fluctuations in the radiation temperature, the grav- 
itational potentials, and the peculiar velocity of the radiation where they last scattered, 
as the temperature anisotropics that we observe today. Only a small fraction of CMB 
photons - current results from CMB polarization measurements [9] indicate around 10 
per cent - underwent further scattering once the universe reionized, due to the ionizing 
radiation from the first stars. 

The details of hydrogen recombination are complicated since it occurs out of quasi- 
equilibrium, but the broad details are well understood [33, 34]. Most hydrogen atoms 
recombined through the (forbidden) two-photon transition between the 2s and 1^ states. 
Direct recombination to Is, or via allowed transitions from excited states, would have 
produced a large excess of photons more energetic than Ly-a (i.e. the 2p-ls transition) 
over those in the Wien tail of an equilibrium blackbody spectrum, and would therefore 
have been rapidly re-absorbed. 

Understanding hydrogen recombination at the level required for sub-percent accuracy 
in CMB observables is both a challenging and ongoing activity. A major step forward 
was made in [35] where a number of approximations made in the original calculations 
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2.2. Thermal history 



were improved upon. However, this improved calculation still misses a number of effects 
that are thought to influence the CMB power spectra by more than 1% for / > 1500 and 
so there has been a flurry of recent activity in this area driven by the imminent launch of 
the Planck satellite. For a recent review of some of the issues involved, see [36]. 



2.3. Linear anisotropy generation: gravity 

We shall mostly be concerned with linearised scalar (density) perturbations about 
spatially-flat models in general relativity. It is convenient to work in the conformal 
Newtonian gauge in which the spacetime metric is 

ds^^a^{n)[{l+2\if)dri^-{l-2^)dx^] , (9) 

where rj is conformal time (related to proper time t by dt = adrj), a is the scale 
factor in the background model and x is comoving position. The physics of the CMB 
is particularly intuitive in this gauge. The two scalar potentials and y represent the 
scalar perturbation to the metric, and play a similar role to the Newtonian gravitational 
potential. In the absence of anisotropic stress in matter and fields, ^ = Xj/. 

We parameterise a CMB photon by its energy e/a as seen by an observer at rest in 
our coordinate system, where e is the comoving energy, and the direction cosines e = 
e^), with = 1, relative to an orthonormal triad of spatial vectors, fl~^(l + 0)<9;. 
The comoving energy is constant in the absence of metric perturbations. The photon 
four-momentum is then 

/7^=a-2£[l-v/,(l+(^)e] . (10) 

We write the photon distribution function, f{x^,p^), as the sum of a background part 
/(e), which is a Planck function with E/f{r\) = e/{af) and so has no explicit time 
dependence ?Lsaf{r\)— const., and a linear perturbation. Denoting small changes in the 
temperature of the radiation for different directions by 7'(77)0(e), the perturbation to the 
distribution function can be written as — 0(e)ed//de so that 

/(77,x,e,e) =/(£)[! -0(T],x,e)dln/(e)/dlne] . (11) 

The quadrupole and higher multipoles of the fractional temperature fluctuation are 

gauge-invariant. 

The distribution function is conserved along the photon path in phase space when 
there is no scattering. Parameterising the path with r\, we have dx/dr\ — yA)e 
from = d;c"/dA (with A an affine parameter). The evolution of e and e follows from 



the geodesic equation; at linear order 

de/dr] = -edyz/dT] -^e(^ + 1//-) , (12) 

de/d7] = -(V-ee-V)((^ + v/) , (13) 

where overdots denote differentiation with respect to r\ and di/z/dr] = i/A + e ■ Vi// is 



the derivative of \\r along the photon path. The comoving energy thus evolves due to 
variation in the gravitational potentials experienced by the photon, and the direction 
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is deflected by the gradient of </) + i/A perpendicular to e (i.e. gravitational lensing). The 
derivative of the distribution function along the spacetime path affixed e and e is sourced 
by the gravitational redshifting in Eq. (12), so that 



d0(77,x,e) 
dr] 



grav 



grav 




(14) 



The gravitational lensing effect does not contribute at linear order since the zero-order 
distribution function / is isotropic. If our observation event has coordinates (t]j?,x/;), 
a photon would have last scattered around recombination at the event (t/*,x*), where 
V* = Vr~ X* ^nd = X/? — ex* ■ On integrating Eq. (14) between these events, the first 
term on the right produces a temperature perturbation given by the difference between 
X]/ at last scattering and at reception. A photon suffers a net redshift if it climbs out of a 
deeper potential well than it is received in and this corresponds to a negative temperature 
fluctuation. The second term on the right of Eq. (14) gives the integrated Sachs-Wolfe 
contribution to 0, 



where the integral is taken along the line of sight. At linear order, the only contributions 
to this integral are from the evolution of the potentials around recombination, due to the 
duration of the matter-radiation transition (the early ISW effect), and the late-time ISW 
effect, when the potentials decay as dark energy becomes dynamically important. For 
the latter, in traversing a potential well that is getting shallower in time (i.e. -|- > 0), 
photons receive a net blueshift and a positive contribution to 0. 



The dominant scattering mechanism to affect CMB anisotropics is classical Thom- 
son scattering off non-relativistic, free electrons (with proper number density Re)- The 
evolution of the temperature fluctuation due to scattering takes the form 



where Or is the Thomson cross-section. The first term on the right describes scattering 
out of the beam; the second, scattering into the beam and involves the familiar 1 -|- cos^ 6 
dependence where is the scattering angle; the final term is a Doppler effect that 
arises from the peculiar (bulk) velocity of the electrons, Vb = dx/dr/. For the latter, 
scattering from a moving electron boosts the energy of a photon scattered in the direction 
that the electron is moving. Note that the effect of scattering vanishes for an isotropic 
distribution function and Vb = 0. If we neglect the anisotropic nature of Thomson 
scattering, replacing 1 -|- cos^ by its angular average 4/3, Eq. (16) reduces to 




(15) 



2.4. Linear anisotropy generation: scattering 



d0 
d77 



scattering 




d0 
drj 



scattering 



aneOT(0 - 00 - e • Vb) , 



(17) 
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where ©o = / Ode/ (4;r) is the monopole of the temperature fluctuation (which equals 
5y/4 with 5y the photon energy density contrast). As expected, scattering tends to 
isotropise the photons in the baryon rest frame, ^ 0o + e ■ Vb- 

Combining the gravitational and scattering contributions gives the Boltzmann equa- 
tion for 0, which is most conveniently written as 

^[e-^(0 + V/)] = -te-^(^v^+e-Vb + Y^ydm0(e,m) [l + (e -111)2]^ 

+ e-^(0 + VA), (18) 

where t(t]) = /^^a/ieOrdT]' is the optical depth to Thomson scattering. The quantity 
— te~^ is the (conformal) visibility function, introduced in Sec. 2.2. 

A formal solution to Eq. (18) can be found by integrating along the line of sight from 
some early time, when T — oo, to the reception event R where T = 0: 

[0(e) + v/]^= TsdT], (19) 





with S the source term on the right-hand side of Eq. (18). The standard Boltzmann 
codes, CMBFAST [37] and CAMB [38] evaluate this integral directly and efficiently 
by expanding the source term in Fourier modes. For our purposes here, we gain useful 
insight into the physics of anisotropy formation by approximating last scattering as 
instantaneous, in which case the visibility function is a delta-function, and ignoring the 
effect of reionization and the anisotropic nature of Thomson scattering; Eq. (19) then 
approximates to 

fR . 

[0(e) + V/]i; = 0o|* + V^* + e-Vb|*+y {^ + ^if)d^. (20) 

The various terms have a simple physical interpretation. The temperature received along 
e is the isotropic temperature of the CMB at the last scattering event on the line of 
sight, 00 1*, corrected for the gravitational redshift due to the difference in the potential 
V'* — Wr^ ^iid the Doppler shift e • Vb|* resulting from scattering off moving electrons. 
Finally, there is the integrated Sachs- Wolfe contribution from evolution of the potentials 
along the line of sight. 



2.4.1. Spatial-to-angular projection 

The main contribution to the CMB temperature anisotropics comes from the projec- 
tion of 00 -I- at the time of last scattering onto the sphere of radius centred on xr. 
For a general, statistically homogeneous and isotropic field /(x), with a Fourier repre- 
sentation 
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the spherical multipoles of the projection of / on the sphere |x — — are 

^^/(k)7K^:^*)iz;(k)e*"^ • (22) 

Here, ji{kx*) are spherical Bessel functions. For large / they peak sharply around kx* = I 
so that the multipoles //„ receive their majority contribution from Fourier modes of 
wavenumber k ~ l/X*- For lii^c of sight, such wavevectors aligned perpendicu- 
lar to the line of sight give angular structure on the sky with characteristic period 
2n I {kx* ) ~ 27r// matching that of F/^, (n) . The ji {kx* ) have oscillatory tails for kx* ^ / 
so that modes with k> l/x* also contribute to the this corresponds to shorter wave- 
length Fourier modes not oriented perpendicular to the line of sight. Introducing the 
dimensionless power spectrum, ^f{k), of /(x), 

(/(k)r(k')) = ^^/(^)5(k-k') , (23) 
the angular power spectrum of the projection is 

c{ = 4kJ ^f{k)jf{kx,)dlnk . (24) 

The contribution to Cf from 0o + at last scattering takes this form with / = ©o + i/a. 

The arguments for the projection of a field that is locally a dipole, like e ■ Vb, are a little 
different. For scalar perturbations, the Fourier components of Vb take the form ikvb(k) 
and so the projection of e • Vb has multipoles 

(e • Vb)/. = -4m^ J j0j-,vMfi{kX*)YL{^)^''-''' ■ (25) 

Here, j'i{kx*) is the derivative of ji{kx*) with respect to its argument. It does not have 
the same prominent peak as ji{kx*), rather having positive and negative wings either 
side of kx* = I. The projection from k to I is therefore much less sharp for terms that 
are locally a dipole compared to monopole terms. Physically, this arises because there 
is no contribution to the projection of e • Vb from Fourier modes with their wavevector 
perpendicular to the line of sight. 



2.5. Acoustic physics 

The photon isotropic temperature perturbation ©o = 5y/4 and baryon velocity Vb 
depend on the acoustic physics of the pre-recombination plasma. Before recombination, 
the comoving mean-free path of CMB photons is 

lp = {ancaT)-^ = -t~^ ^4.9 xlO^^ia^h^y^il+zy^Mpc . (26) 

On scales larger than this, Thomson scattering keeps the photons isotropic in the baryon 
rest frame and the photon-baryon system can be modelled as a tightly-coupled fluid 
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characterised by a single peculiar velocity Vb- The Euler equation for this system, in the 
Newtonian gauge, is 



+ TT^-^^^ + 4(TT^^^^+ ^'^ = ^ ' ^^^^ 

where R = Pb/{Py + Py) = 3pb/ (4py) °^ a and p,- and are the unperturbed energy den- 
sity and pressure of species i. The drag term, involving the conformal Hubble parameter 
J^f = d/a, arises since the peculiar velocity of non-relativistic massive particles (the 
baryons) decays as l/a in an expanding universe. The V5y term comes from the gradi- 
ent of the photon pressure. This acts on the full inertia of the photon-baryon system and 
so its effect reduces for increasing R. The final term describes (universal) gravitational 
infall in th potential gradient Vi/A^. The energy density of the photons is not changed by 
scattering off cold electrons and so 5y evolves as for a non-interacting photon gas: 

5r+^V-Vb-40 = O. (28) 

The second and third terms in this equation are the usual p+p times volume expansion 
rate, where the perturbed expansion rate receives contributions from V • Vb and the 
perturbation — 3^d^x to the volume element (see Eq. (9)). 

Combining Eqs (27) and (28) gives the equation of a damped, harmonic oscillator 
driven by gravity [39]: 

g ^ v2gy = 4(i)- + ^^(^ + -VV- (29) 

^ l+R ^ 3{l+R) ^ ^ l+R^ 3 ^ 

First, consider the solutions to the homogeneous equation, ignoring the effects of gravity. 
On sub-Hubble scales, k J^, the WKB approximation can be used to show that 

(l+i?)-V4cosA:rs and 5y^ {l+Ry^/'^sinkr, , (30) 



are independent solutions, where the sound horizon = Jq drj' / ^/3{l+R). These are 



acoustic oscillations at instantaneous frequency k/ ^y3{l +R), where l/y3(l+7?) is 
the sound speed of the photon-baryon system, with a slow decay in amplitude. We shall 
show shortly that for adiabatic initial conditions, as expected in single-field models of 
inflation, the system oscillates essentially in the cosine mode. All Fourier modes of a 
given k, irrespective of the direction of their wavevector, reach extrema of their acoustic 
oscillations at the same time. This phase coherence is an inevitable consequence of 
the dominance of a single "growing mode" in the primordial perturbation at the end of 
inflation, and is not present in models that continually source perturbations such as those 
with topological defects [40] . Phase coherence implies a series of acoustic peaks in the 
CMB angular power spectrum (see Fig. 2) at multipoles / such that l/x* corresponds 
to the wavenumber of modes that have reached extrema at i]*. For adiabatic initial 
conditions, the acoustic peaks are located approximately at / = nKX^/r^if]*). 
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2.5.1. Initial conditions and gravitational driving 



We can get the initial conditions for the oscillator deep in the radiation era by con- 
sidering the self-gravitating perturbations of a radiation fluid. Such a treatment makes a 
number of assumptions. First, neutrinos are treated as a fluid but, after decoupling, they 
actually free-stream and so do not cluster appreciably on sub-Hubble scales. Moreover, 
free-streaming generates anisotropic stress ~ pv{kr})^^ in the neutrinos which produces 
a time-independent 14% offset between (the constant) and i/a on large scales. Second, 
we are assuming that the perturbations in the baryonic matter and CDM have a negligi- 
ble effect on the gravitational potentials. This is appropriate for adiabatic initial condi- 
tions deep in the radiation era, where all species are initially distributed homogeneously 
on large scales on the same spacetime hypersurfaces, i.e. the density perturbations of 
all species vanish on hypersurfaces over which the total density perturbation vanishes. 
In any other gauge, P/5;/p; is equal for all species and there are therefore no relative 
perturbations in the number density of particles. For adiabatic initial conditions, the 
total energy density contrast in the radiation era is dominated by the radiation perturba- 
tions. This should be compared to isocurvature initial conditions, where, for example, 
the CDM density contrast is initially much larger than the radiation contrast in such a 
way as to give no net curvature perturbation. 

In a radiation fluid, Jif = l/r] and the trace of the ij Einstein equation (see, for 
example, [41]) gives 

i^ + -(j>--L(j>^-l-dr, (31) 

where the subscript r denotes radiation (photons plus neutrinos). The 00 Einstein equa- 
tion is 

V'0 = ^(^0 + 0) + ^5,, (32) 

so the potential evolves as 

+ ^0-^v20 = O, (33) 

with regular (growing-mode) solution ^ <x ji{kr} /\^) /{krj). This is constant outside 
the sound horizon, krj <^ \/3, but undergoes decaying oscillations on small scales with 
asymptotic behaviour^ ^ <x —\/3cos{kr] / y/S) / {kr])^. Solving for 5r from Eq. (32) gives 
a solution with 5r = —20(0) on large scales (i.e. 5r is constant) and constant amplitude 
oscillations, 5r = — 6cos(A:r//-\/3)0(O), well inside the sound horizon. As advertised 
above, adiabatic initial conditions excite the cosine mode of the oscillator. Note also that 
the radiation starts off over-dense in potential wells. 

If we now return to Eq. (29), we can make two important observations that affect the 
height of the acoustic peaks in the CMB power spectrum. First, for modes that enter 
the sound horizon between matter-radiation equality and recombination, the oscillations 
occur in potentials that approximately equal the primordial value 0(0). More carefully. 



The presence of the CDM, which does not undergo acoustic oscillations, leads to a A:-dependent offset, 
~ k~^lnk, in the oscillation as the matter-radiation transition is approached. 
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FIGURE 3. Contribution of the various terms in Eq. (20) to the temperature-anisotropy power spectrum 
from adiabatic initial conditions: 5y/A + (denoted SW for Sachs-Wolfe [42]; magenta); Doppler effect 
from Vb (blue); and the integrated Sachs-Wolfe effect (ISW; green) coming from evolution of the potential 
along the line of sight. The units of the spectrum are arbitrary. 



the potential decays by 10% through the matter-radiation transition on large scales - 
although this is not complete at recombination - and, by Eq. (28), dy is enhanced by 
20%. The potential shifts the midpoint of the oscillation to 5j ~ —4(1 +R)^f which 
corresponds to the over-density needed for photon pressure to balance gravitational 
infall. Since baryons contribute to the inertia but not the pressure, they enhance 8y at 
the midpoint. The source term for the temperature anisotropy, 0o + = 5y/4 + ^f, 
therefore oscillates about — with an amplitude of 30(0)/ 10 (ignoring non-zero R 
in the amplitude). The dependence of the midpoint on the baryon density makes the 
relative heights of the first few acoustic peaks sensitive to the baryon density. The 
second observation is that for modes shorter than the sound horizon at matter-radiation 
equality, the decay of the potential during the oscillations in the radiation renders the 
driving term on the right of Eq. (29) negligible during the matter era. However, our 
treatment of the radiation fluid above shows that the driving term acts resonantly around 
the time of sound-horizon crossing in the radiation era and enhances the amplitude of the 
acoustic oscillation of 5y/4 towards the asymptotic value 30(O)/2, five times the value 
on large scales. Increasing the matter density limits the scales for which resonant driving 
is effective to smaller scales. In practice, diffusion damping (see Sec. 2.7) exponentially 
damps the amplitude of the oscillations on small scales and the asymptotic amplitude is 
not attainable. 
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2.6. Cj for adiabatic initial conditions 

The contributions of the various terms in Eq. (20) to the temperature anisotropy power 
spectrum for adiabatic initial conditions are shown in Fig. 3. The main contribution is 
from 00 + y/. The plateau on large scales comes from modes that have not had time to 
oscillate by recombination; for such modes, 0o + 30(O)/1O = 0/3. It is convenient 
to express the initial value of the potential in terms of the comoving-gauge curvature 
perturbation, 

M=-6-^^^±^ (34) 
^ 47tGa^{p+p)' ^^^^ 

which is conserved on super-Hubble scales for adiabatic perturbations independent of 
the equation of state relating the total pressure p and energy density p. During radiation 
domination, =^ = —30(0) /2 ignoring neutrino anisotropic stress, and so on large scales 
@0+Y^ —^/5. This directly relates Cf to the primordial power spectrum ^^{k) 
(which is calculable given a model of inflation) on large scales. In particular, a nearly 
scale-invariant ^^{k) projects to give /( / + 1)C/ const. The prominent acoustic peaks 
in Fig. 3 reflect the acoustic oscillations of the photon-baryon system. The first peak 
corresponds to modes that have just reached the first maximum of their oscillation by 
last scattering. The resonant-driving arguments given earlier suggest that the acoustic 
peaks should grow in amplitude with increasing /. In practice, an exponential decline in 
power is seen because of diffusion damping. 

The late-time ISW effect adds power incoherently to the contribution from recombi- 
nation since the former is sourced at late times. However, the early-ISW effect, which 
peaks close to the first acoustic peak, adds power coherently [39]. It enhances the first 
peak considerably since <]) + yf has the same sign, i.e. positive in potential wells, as 
00 + V- It is Isss important for the higher peaks since the potential at last scattering for 
these is suppressed by the acoustic oscillations during radiation domination. 

The baryon peculiar velocity oscillates 7t/2 out of phase with 5y by Eq. (28), and 
so the Doppler term tends to fill in the troughs of the angular power spectrum. The 
oscillatory structure would be destroyed completely were it not for two effects: the 
geometric effect noted for the projection of dipole-like fields in Sec. 2.4, and the effect 
of baryons. The latter shift the midpoint of the oscillations m@Q + \]/ from zero, so its 
zeroes do not coincide with the extrema of the peculiar velocity. Furthermore, baryons 
reduce the amplitude of the oscillations in Vb by a factor of (1 +R). 



2.7. Complications to the simple picture 

2.7.1. Photon diffusion 

On small scales, the CMB anisotropics are exponentially damped due to photon 
diffusion. Before recombination, fluctuations with wavelengths comparable to or smaller 
than the mean-free path of photons to Thomson scattering are damped out, as photons 
can diffuse out of over-dense regions into neighbouring under-densities. How far will a 
photon have diffused by last scattering? To answer this, consider the comoving mean- 
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free path of the photons, Ip, which is related to the Thomson cross section and the 
number density of free electrons via 



lp = -^. (35) 

In some interval of conformal time, dr\, a photon undergoes drj/lp scatterings and 
random walks through a mean-squared distance Ipdr] /Ip — Ipdr] . The total mean-squared 
distance that a photon will have moved by such a random walk by the time T]* is therefore 

(36) 

atieOj kj) 



which defines a damping scale kj^^. The photon density and velocity perturbations at 

wavenumber k are damped exponentially by photon diffusion, going like e^/^o; this 
produces a similar damping in the angular power spectrum Q (see Fig. 3). Evaluating 
kn around recombination gives a comoving damping scale ~ 30Mpc. 



2.7.2. Reionization 

Once structure formation had proceeded to produce the first sources of ionizing 
radiation, neutral hydrogen and helium began to reionize. The resulting free electron 
density could then re- scatter CMB photons, and this affected the observed CMB in 
several ways. The only significant effect for the temperature anisotropics in linear theory 
is a uniform screening by q^'^'^, where Tre is the optical depth through reionization, 
on scales / > 10. Anisotropies on such scales are generated from perturbations that 
are sub-Hubble at the time of reionization. For these modes, the radiation is already 
significantly anisotropic at reionization since the perturbations vary significantly over 
the scattering electron's own last-scattering surface. The net effect of in-scattering from 
different lines of sight thus averages to zero, leaving a suppression by e~^ due to 
scattering out of the line of sight. In the power spectrum, this becomes e~^^"=. However, 
Fourier modes that are still super-Hubble at reionization produce negligible anisotropy 
by the time of reionization since the perturbations over the scattering electron's own last 
scattering surface are almost uniform. Scattering isotropic radiation has no effect, and so 
reionization does not alter the temperature anisotropies from recombination for / < 10. 
Scattering around reionization also generates new large-angle polarization - see Sec. 3.3 
- and, at second-order, small-scale temperature anisotropies; see Sec. 6.2. 

Note that screening at reionization makes the determination of the amplitude and 
shape of the primordial power spectrum from temperature data alone rather degenerate 
with Tre- 
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2.73. Gravitational waves 



Tensor modes, describing gravitational waves, represent the transverse trace-free per- 
turbations to the spatial metric: 

d*^ = fl2 ) [dT]2 - + hij)dJdoi^\ , (37) 

with h\ — and = 0. A convenient parameterisation of the photon four-momentum 
in this case is 



2 ^ 



(38) 



where, as for scalar perturbations, = 1 and e is a times the energy of the photon as seen 
by an observer at constant x. The components of e are the direction cosines of the photon 
direction for this observer on an orthonormal spatial triad of vectors aT^idi — h\dj/2), 
and are constant in the unperturbed universe. The evolution of the comoving energy, e, 
follows from the geodesic equation: 

~ + 5%-V=0, (39) 

and so the Boltzmann equation for the fractional temperature anisotropics is 

^+e-V0 = -anegxQ+ ^T/^'^ / dih0(£,ih) ri + fe-ih)^] 
dr\ \6k J / L V /J 

-Uij^eL (40) 



All perturbed scalar and vector quantities vanish for tensor perturbations so 0(e) has 
only / > 2 moments. Neglecting the anisotropic nature of Thomson scattering, or, equiv- 
alently, the temperature quadrupole at last scattering, the solution of Eq. (40) is an inte- 
gral along the unperturbed line of sight: 

[®{'e^)]R=-\j\-%je'eUr] . (41) 

The physics behind this solution is as follows. The time derivative hij is the shear 
induced by the gravitational waves. This quadrupole perturbation to the expansion rate 
produces an anisotropic redshifting of the CMB photons and an associated temperature 
anisotropy. 

The shear source term for the tensor anisotropics is locally a quadrupole since hij is 
trace-free. Projecting this contribution at a given time onto angular multipoles on the sky 
is a little more involved than for scalar perturbations [43]. We start by decomposing hij 
into circularly-polarized Fourier modes, 

% = E/pS72CWe'''^ (k) = -Lmf >(k)ft(^^)(k) . (42) 
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("±2") 

Here, the basis tensors m-^ '^(k) are complex, symmetric trace-free and orthogonal 
to k, and /i^^^) (k) are scalar Fourier amplitudes. For k along the z-axis, we choose 
m-^^^ = (x ± jy),(x ± iy)j/2. For such Fourier modes the shear source is 



Aif ) (^z). V - ^ [(X ± iy) • e]' = ^ sin2 ee^^'^ = ^ ^y2±2(e) , 
The projection of this at comoving distance % is 



(43) 



^'^''y2±2(e)e-'^^-^« = 4^£(-/)^V2LTT;L(^;t:)y2±2(e)yLo(e) 



{-i)\lL+\)jL{kx) 




X 



,(44) 



where we have used the Rayleigh plane-wave expansion. The 3 j symbols arise from 
coupling the angular dependence of the source to that of the plane wave. Since they 
force L and / to have the same parity, the /th multipoles of the projection involve the / 
and / ± 2 multipoles of the plane wave. We can simplify further by writing out the 3 j 
symbols explicitly and using the recursion relations for Bessel functions to express ji^2 
in terms of 7/. The final result is remarkably compact: 



' {l + 2)\ji{kx) , 
(/-2)! 



'Z±2 



(e). 



(45) 

The spatial-to-angular projection is now controlled by j i{kx) / (kx)^ ■ As for scalar 
perturbations, this is concentrated on multipoles I ^kX- The result for general k can be 
obtained by a rotation Ok,0), where 9^ and 0k are the polar and azimuthal angles 
of k and the third Euler angle can be taken to be zero. Pre-empting a little the discussion 
in Sec. 3, we can relate the rotation matrices to spin-weighted spherical harmonics, 
(0k,0k,O) 



^m±2 



^4K/{2l+l)^2Y;jk), SO that 



(k)eVe-'^'^-^ = -V2^A(±2)(k) £(-/)' 

Im 



{1 + 2)1 jijkx) 

(/-2)! (kx)^^'^i'nimm{e). (46) 



This generalises the scalar result for which the multipoles were given in Eq. (22). 

Figure 4 compares the power spectra due to gravitational waves with those from 
scalar perturbations. The former is shown at the maximum amplitude allowed by current 
data. Note that the spectra from gravitational waves fall sharply on scales that are sub- 
Hubble at recombination (/ > 60) since the amplitude of gravitational waves decays 
away as 1 /a inside the Hubble radius. This limits the power of temperature anisotropies 
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FIGURE 4. Temperature (black), is-mode (green), B-mode (blue) and T-E cross-correlation (red) 
CMB power spectra from scalar perturbations (left) and tensor perturbations (gravitational waves; right). 
The amplitude of the tensor perturbations is shown at the maximum amplitude allowed by current data 
(r = 0.22 [44]). The B-mode spectrum induced by weak gravitational lensing is also shown in the left-hand 
panel (blue; see Sec. 6.1.2). 



to constrain gravitational waves since the sampling variance of the dominant scalar 
perturbations is large at low /. Fortunately, CMB polarization provides an alternative 
route to detecting the effect of gravitational waves on the CMB which is not limited by 
cosmic variance [45, 46]; see Sec. 3. 



2. 7.4. Isocurvature modes 

Adiabatic fluctuations are a generic prediction of single-field inflation models. How- 
ever, multiple scalar fields typically arise in models inspired by high-energy physics, 
such as the axion model [47], curvaton [48] and multi-field inflation [49, 50]. In such 
models, if the fields decay asymmetrically and the decay products are unable to reach 
chemical equilibrium with each other, an isocurvature contribution to the primordial 
perturbation will result. The simplest, and best-motivated, possibility is an isocurva- 
ture mode where initially the dominant fractional over-density is in the CDM, with a 
compensating (very small) fractional fluctuation in the radiation and baryons [51]. The 
amplitude of the CDM isocurvature mode is quantified by the gauge-invariant quantity 
^ = 5c — 35y/4, where 5c is the CDM fractional over-density. Generally, can be cor- 
related with the curvature perturbation 3$, for example in the curvaton and multi-field 
models. 

In the CDM mode, the photons are initially unperturbed, as is the geometry: 5y(0) = 
= ^(0) and Vb = 0. The different equations of state of the CDM and radiation lead to 
the generation of a curvature perturbation. On large scales, ^ grows like a in radiation 
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domination and is constant in matter domination with an asymptotic value y /3. 
The gravitational potential follows suit and, for super-Hubble modes at last scattering, 
0* = —3^%-^/5 = —y/5. On large scales dy/A = ^ is preserved, so that the source term 
for the temperature anisotropics on large scales is 0o + ~ —2^/5. For a general mix 
of isocurvature and adiabatic perturbations, the large-angle temperature anisotropics are 
related to the primordial perturbations by /ST /T ~ —M/5 — 15^ / 5. 

The evolution of the potential for isocurvature modes makes the driving term in 
Eq. (29) mimic the sine solution of the homogeneous equation, and so dy oscillates 
as sin^rs inside the sound horizon. This n/2 phase difference compared to adiabatic 
initial conditions shifts the acoustic peaks to multipoles / = {n + \/l)Ttx*/'''s{'^*)- The 
different peak positions for isocurvature initial conditions allow the CMB to constrain 
their relative contribution to the total fluctuations. The constraints on isocurvature modes 
depend strongly on the assumptions made about the shape of the primordial power 
spectrum for and any correlations with M. For uncorrelated, axion-type models with 
scale-free CDM isocurvature fluctuations, the ratio of the primordial spectra for to 
M <0.2 from CMB data alone [44]. Recent constraints in more general models can be 
found in [52]. 

3. INTRODUCTION TO CMB POLARIZATION 

As the universe recombines, the mean-free path to Thomson scattering grows and tem- 
perature anisotropics start to develop. Subsequent scattering of the CMB quadrupole 
generates partial linear polarization in the CMB [5]. The signal is small, with r.m.s. 
~ 5 juK, but contains complementary information to the temperature anisotropics and is 
therefore being actively pursued by observers. 



3.1. Polarization observables 

Polarized specific brightness is most conveniently described in terms of Stokes pa- 
rameters /(n; v), 2(n; v), U (n; v) and V(n; v). The / parameter is the total intensity, Q 
measures the difference in brightness between two orthogonal linear polarizations, and 
V measures the circular polarization. The U parameter, like Q, measures linear polar- 
ization and is defined as the difference in brightness between two linear polarizations at 
45° to those used to define Q. Thomson scattering of the CMB does not generate cir- 
cular polarization so we expect the primordial V — Q. The Stokes parameters Q and U 
are functions of frequency as well as direction. In linear theory, their frequency depen- 
dence follows that of the temperature anisotropics; dividing by the derivative of the CMB 
blackbody (i.e. Planck) brightness with respect to (log) temperature allows us to express 
the Stokes parameters as frequency-independent fractional thermodynamic equivalent 
temperatures. We shall denote the equivalent temperatures by Q{n) and U (n). 

The linear Stokes parameters depend on the choice of basis. For a line of sight n, we 
define the Stokes parameters on a local x-y basis defined by —0,(j) of the spherical-polar 
basis at n. Note that — 0,<^ and — n = e - the radiation propagation direction - form 
a right-handed basis. The reader is cautioned that this convention differs from some 
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of the CMB literature, most notably [10]. Under the right-handed rotation about the 
propagation direction, ± /0 — > e '^(0 ± the complex polarization 

Q±iU ^ t^^'^^iQ ± iU) , (47) 

and so 2 ± iU is said to be spin =f2. The complex polarization can be expanded in terms 
of a real scalar field and a real pseudo-scalar field L/m^/mi'/mlA)' where 

the summations are over / > 2, as [10] 

(2 ± iU) (n) = Y,{Eim T iBim)T2Yim{n) . (48) 

Here, the ^2Jlm(n) are spin-^2 spherical harmonics [53]. Spin-5 spherical harmonics 
are defined by 





1^1)!. 







where S"'*' = (— l)*8l*L The spin-raising and lowering operators, 5 and 6 have action 

g^T] = -sin*0((90-|-/cosec(9^)(sin~*0iT]) (50) 
g^T] = - sin~* 6 {de - /coseci?^ ) (sin* 0,7]) (51) 

when acting on a spin-* quantity ^77. They raise and lower the spin by one unit respec- 
tively. 

The decomposition of the polarization into E and B modes is analogous to writing a 
vector on the 2- sphere, Va, as the sum of a gradient part, VqVe, and a divergence-free 
part, which can always be written as Ca^V/^Vg where Eat is the alternating tensor. The 
components of Va on the null basis m± = (0 ± z^)/ \/2 are then 

V+ =\.m+ = j^{d0+ icosccd)VE -^{de + /cosec0)ys = -^^{Ve - iVe) , 
y_ =\.in_ = ^{d0- /cosec0)y£ + ^ (^0 - Jcosec0)ys = -j^^Vs + iVe) (52) 

We see that the spin-1 V+ is written as the action of the spin-raising operator on the 
potential —Ve + /Vg whose real and imaginary parts are the (real) E and B potentials 
for V. Things work similarly for polarization since Q ± iU are the null components of a 
symmetric, trace-free tensor [54]. 

To gain further intuition for E modes and B modes, consider the case where the 
potentials Y^imEimYimi^) Slid Y^imBimYimi^) behave locally like a plane wave across a 
small patch of the sky (that can accurately be treated as flat). The E and B contributions 
to the polarization are depicted in Fig. 5. Quite generally, a given potential i// generates 
a B mode with Q^iU ^ -id^ Xj/ and Q-iU = id^ Xj/, and an E mode with Q + iU ^Q^xi/ 
and Q — iU = d^ij/. Since / = e^'^/^, the 5-mode polarization pattern is obtained from 
the E mode by rotating by 45° about the line of sight. 

Statistics of CMB polarization. The decomposition of the polarization field into E 
and B parts is invariant under rotations, and the E and B multipoles transform like 
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FIGURE 5. Polarization patterns for an E mode (left) and B mode (right) on a small patch of the sky 
for potentials that are locally Fourier modes. The shading denotes the amplitude of the potential and the 
headless vectors denote the unsigned direction of the polarization. For the electric pattern the polarization 
is aligned with or perpendicular to the Fourier wavevector depending on the sign of the potential; for the 
magnetic pattern the polarization is at 45 degrees. 



those for the temperature anisotropies, e.g. Jim' ^mm'^im'- Under the operation 

of parity, {Q±iU){n) — > (2=f/[/)(— n) so that (— l)'^/^ (electric parity) while 

Bim (magnetic parity). 

Rotational and parity invariance in the mean limit the non-zero two-point correlations 
between polarization (and temperature) multipoles to 

{EimEp„^/) = Siifd„,„/Cf (53) 
{BimB}^/) = Siifd„,„/Cf (54) 
{TimEp^f) = Siifd„,„/Cj^ , (55) 

which define the power spectra Cf , Cf and Cj^ . Note that there is no correlation 
between B and either AT or E. 

Assuming Gaussian statistics, the cosmic variance in the polarization (auto-)power 
spectra are 

-^iCf) = :^^iCf)\ var(Cf) = ^(Cf)^ (56) 
For the cross-spectrum, the cosmic variance is only slightly more complicated [45, 46]: 

var(Cf^) = [{Cn^ + Cfcf] . (57) 

The second term arises from chance correlations between AT and E in the single 
realisation of the sky that we have available. It is the only term present in the cosmic 
variance of estimators for the cross-spectra that vanish in the mean by parity, Cf^ 
and Cf^. Due to the presence of correlations between the temperature and the electric 
polarization, estimates of Cf and Cf have non-vanishing covariance, and each is also 
correlated with the estimator for Cf^; see [45, 46] for details. 

The form of the higher-order polarization spectra implied by rotational invariance in 
non-Gaussian theories is similar to that for temperature and is described in [55]. 
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3.2. Thomson scattering and free-streaming 



Polarization is generated by scattering the quadrupole of the temperature anisotropy 
around recombination and reionization. Prior to recombination, on scales larger than 
the mean-free path, Thomson scattering keeps the CMB radiation very nearly isotropic 
in the rest-frame of the baryons, and therefore unpolarized. It is only as the mean-free 
path grows as recombination is approached that the radiation can free-stream over an 
appreciable fraction of a wavelength of the perturbations between scatterings and so 
generate a temperature quadrupole. Subsequent scattering of this quadrupole generates 
linear polarization. After the CMB last scatters around recombination, the polarization 
is preserved until further scattering at reionization (see Sec. 3.3). 

Thomson scattering changes the linear polarization state of the radiation as [43]^ 

d(!2±/C/) = td77(!2±/f/)(e)-^td77 V ( E2m - -^&2m) ilYlmi^) (58) 

5 \m\<2\ V6 J 

Here, 0/„ are the multipoles of the fractional temperature fluctuation, 0(7], x,e) = 
£;„0;m(T7,x)y;„(e), and, recall, t = — aneCTx. The first term arises from scattering out 
of the beam and reduces the polarization. The second term, arising from in-scattering, 
involves the quadrupoles of the temperature anisotropics and the E-mode polarization. In 
particular, scattering of unpolarized radiation generates linear polarization d{Q + iU) = 
-tdri 52 @2mY2m) /20 which is an E-mode quadrupole. However, the polarization 
does not generally remain in this state in the presence of spatial inhomogeneities in 
the polarization induced over the last scattering surface. To see how this works, it is 
simplest to consider the important cases of density perturbations and gravitational waves 
separately. 



3.2.1. Density Perturbations 

Density perturbations transform as scalars under reparameterisation of spatial coordi- 
nates. In Fourier space, this means that every perturbed tensor is azimuthally-symmetric 
about the wavevector k. The same is true of the temperature fluctuation and so for the 



When discussing the production and propagation of polarized radiation, it is convenient to associate the 
radiation with the point e on the sphere. At this point, the propagation direction is outwards so a right- 
handed polarization basis is formed by 6 and ^ there; this is equivalent to using 0,— ^atn = — e and so 
the Stokes parameters are the same in both descriptions. However, {Q + iU){e) is then spin +2 and its 
appropriate expansion in spin harmonics is 

{Q ± iU) (e) = ± iBi„)±2Yi„ (e) . 

Im 

The multipoles E/^ and B;^ here are related to those in the line-of-sight description by factors of (— 1)' 
and (— 1)'"*"^ respectively and so the parity-invariant power spectra are unchanged. 
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quadrupole we must have 



©2m(T?,k)F2^(e) ocP2(k-e) = ^ L i^2m(e)F|^(k) , (59) 

\m\<2 \m\<2 

where, in the second equaUty, we have used the addition theorem for spherical harmon- 
ics. As a concrete example, consider scales which are large compared to the mean-free 
path Ip. Applying Eq. (20) over one scattering time, a quadrupole temperature anisotropy 
builds up most efficiently by the Doppler term. If we locate ourselves at the origin, the 
temperature anisotropy is 

0(77 , 0, e) - e • Vb(?] - Ip, -Ipe) fti e • Vb(?7 - /p, 0) - /^e • (e • Vvb) (?] - Ip, 0) . (60) 

The quadrupole part of this comes from the last term which, for a single Fourier mode, 
gives 

£ 02^(T?,k)y2m(e) ~ |/:/pP2(k-e)vb(k) . (61) 

\m\<2 

Note that this is 0{klp) whereas the monopole source terms in Eq. (20) can only 
produce a quadrupole of 0{k^lp). If we choose coordinates with the z-axis along k, 
the temperature anisotropies are all in the m = modes for density perturbations. 

Scattering of the quadrupole generates Q±iU ~ L|m|<2i2*m(k)±2l2m(e) locally. If we 
take k along z, we have Q sin^ d and U = 0. How does the character of the observed 
polarization change by free- streaming (ignoring reionization for the moment)? What we 
observe from this single Fourier mode still has U = 0, but the angular dependence of Q is 
further modulated by the spatial dependence of the perturbation over the last-scattering 
surface. This modulation transfers polarization from / = 2 to higher multipoles (with 
most appreciable power appearing at Z = kx*), but preserves the electric character of the 
polarization [45, 46]. We can see that this is reasonable by noting that the modulated 
polarization field has its polarization direction either parallel or perpendicular to the 
direction in which the polarization amplitude is changing [56]. In more detail, if we take 
k along the z-axis, the observed polarization at the origin is 

iQ±iU)ie) = ((2±/C/)(T?*,-Z*e,e) - (1 -ju2)e-'^^*^, (62) 

where /i = k • e = cos 6. (Note that U is zero in this orientation.) Since the polarization 
is azimuthally symmetric, the same will be true of any E and B modes. Taking 

Q + iU = d^{\j/E + i\lfB), Q-iU^B^{\i/E-i\lfB), (63) 

where y/E = L/>2 •>/(/ -2)1/ (I + 2) !£'/oi/o(0, 0) and similarly for Yb, to satisfy Eq. (62) 
we must have 

{l-li^){x/^±ix/^) oc (1 -ju2)e-*^*'^ , (64) 

where primes denote differentiation with respect to ;U. It follows that i/Zg = (the non- 
zero solutions of = are combinations of / = and / = 1 modes and cannot generate 
polarization) so there are no B modes. Since E and B modes transform irreducibly 
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under rotations, there will be no B modes for a general k. Dropping the proportionality 
constant, the relevant solution of ij/^ = q~^^^*^^ is ao{kx*) + a\ {kx*)fJ- - e~'^^*^ / (kx*)^ 
where ao and a\ should be chosen to remove the / = and / = 1 modes. Expanding the 
exponential gives 



The multipoles for a general direction of k can be found by the appropriate rotation. 
Note that the spatial-to-angular projection of scalar polarization is controlled by the 
same function, ji{kx*)/{kx*)^^ as for the temperature anisotropies from gravitational 
waves (see Eq. (45)). 

The polarization power spectra produced by adiabatic density perturbations are plot- 
ted in the left panel of Fig. 4. They were computed with the Boltzmann code CAMB [38]. 
The E-mode power peaks around / ~ 1000, corresponding to the angle subtended by the 
width of the visibility function at recombination. On larger scales the polarization probes 
the electron-baryon velocity at last scattering, as described above. Acoustic oscillations 
of the plasma prior to last scattering imprint an oscillatory structure on the angular power 
spectra of the CMB observables. Asymptotically, the photon-baryon overdensity oscil- 
lates as a cosine, and the peculiar velocity of the plasma as a sine, to satisfy the continuity 
equation. The peaks in the £'-mode spectrum therefore coincide with the troughs of the 
temperature spectrum. Note also that the temperature and E-mode polarization are cor- 
related [57] and that modes caught at either the midpoint or extrema of their oscillation 
at last scattering give zeroes in the temperature-polarization cross-correlation. Large- 
angle polarization from the last-scattering surface is very small since the generation of 
a quadrupole in the temperature anisotropy occurs via causal physics (photon diffusion) 
that is suppressed outside the horizon. The increase in polarization on large scales in 
Fig. 4 is due to reionization [6] and is discussed further in Sec. 3.3. 



3.2.2. Gravitational Waves 

For gravitational waves, a temperature quadrupole builds up over a scattering time 
due to the shear of the gravitational wave: 



0(e) --^/pVV. (66) 



Consider a circularly-polarized gravitational wave with k along the z-axis so that the 
temperature quadrupole is a |m| = 2 mode. Equations (43) and (58) show that Thomson 
scattering of this temperature quadrupole produce polarization that locally has 

1 1 9,71: 

Q±iU^ Iph^P^ (kz) - J —±2Y2p (e) , (67) 
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where h^P^ (k) is the Fourier amplitude of the gravitational wave with p — ±2 labelling 
the helicity.'* Unlike scalar perturbations, the local polarization has both Q and U non- 
zero and modulating these with the plane wave e"'^*'''^ will produce both E and 5-mode 
polarization [45, 46]. In detail, the observed polarization is 

(<2±/C/)(e) oc ±2F2p(e)e-'^^*'^°^« 

= ±2i2p(e)£v^v^r+T(-/)^Lte)iio(e) 



V2L + TH)'^jLikx.) 



X 



±2^/^(6) 



I L / , 



(68) 



where we have used an addition theorem for the spin spherical harmonics, and the 
functions 



e/(x) 



1 
4 



ftW = 2 



dji{x) 



X dx 



6x +1^^^'^ 



(69) 
(70) 



follow from substituting the explicit form of the 3 j symbols and using the recursion 
relations for the the spherical Bessel functions. Reinstating numerical factors from 
Eq. (67), and rotating to a general k, the non-zero E- and 5-mode multipoles are 



Elm — 
Blm = 



^/9~2 



(71) 
(72) 



We see that gravitational waves project onto both E- and 5-mode polarization. The 
projection from last scattering peaks at multipoles / ~ kx* and divides power roughly 
equally between E and B [43]. Note also that under a parity transformation, A*^^^) (k) — > 
hi^^){-k). Given that ,y/„(-k) = (-l)'_,vy/„(k), Eqs (71) and (72) manifestly give 
Elm {-'^yEim and5/„ -(-l)'5/„, as required, under parity. 

The polarization power spectra from a scale-invariant background of gravitational 
waves with tensor-to-scalar ratio r = 0.22 is shown in the right-hand panel of Fig. 4. 



A more careful treatment of Eq. (40), balancing in- and out-scattering in the presence of the shear source, 
and properly including the polarization dependence of Thomson scattering, introduces a prefactor of 1 0/ 3 
in these expressions; see [43, 58, 59]. 
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As for temperature anisotropies, the spectra peak on scales corresponding to the horizon 
size at recombination. Note that the gravitational wave contribution to the temperature 
anisotropies and electric polarization is sub-dominant to the scalar perturbations, but 
gravitational waves dominate the 5-mode polarization on large angular scales. (On small 
scales, the non-linear 5-modes induced by the action of gravitational lensing on the 
largely E-mode primary polarization are dominant.) The increase in polarization power 
on the largest scales in Fig. 4 is due to reionization, which we now briefly discuss. 

3.3. Polarization from reionization 

For polarization, re-scattering suppresses the signal from last-scattering by a factor 
but also produces a new large-angle signal [6] (see Fig. 4). From Eq. (58) we 
see that the strength of this large-angle polarization varies almost linearly with Tre and 
the size of the temperature quadrupole at reionization. Since reionization appears to 
have completed before dark energy came to dominate the expansion, the temperature 
quadrupole at last scattering is sourced by the Sachs-Wolfe and Doppler effects on 
the last scattering surface of the re-scattering electron. Therefore, the quadrupole is 
dominated by modes with ^(rjre — T]* ) ~ 2. After scattering, the newly-generated l = 2E- 
mode polarization free-streams to give £'-mode (and, for gravitational waves, 5-mode) 
polarization that peaks at multipoles / ~ 2(770 — 77re)/(^re — ^?*)- The position of the 
reionization feature is thus controlled by the epoch of reionization. The detailed history 
of reionization does generate some further small features in large-angle polarization but 
extracting this is rather limited by cosmic variance [60, 61]. 

Currently WMAP is the only polarized dataset with sufficient sky coverage to extract 
Tre. The latest results give Tre = 0.087 ±0.017 [9]. For instantaneous reionization, this 
limits the redshift Zre > 8.2 at 95% confidence, which, combined with quasar absorption 
spectra that indicate reionization was complete by z = 6 (e.g. [62]), points to reionization 
being an extended process. The determination of Tre required aggressive cleaning of 
polarized Galactic foreground emission by the WMAP team, but their result is stable to 
variations in the details of the cleaning process [9]. 



3.4. Current measurements of polarization power spectra 

In Fig. 6, we show all current measurements of the E-mode power spectrum and the 
T-E correlation. WMAP [29] provides the best measurements on large scales while the 
recent QUaD data [63] is most constraining on intermediate and small scales. The data 
is still rather noisy but the qualitative agreement of the polarization spectra to the best- 
fit ACDM model (which is still driven by the higher signal-to-noise temperature data) 
is striking. This is an important test of the structure formation model: the polarization 
mainly reflects the plasma bulk velocities around recombination and these are consistent, 
via the continuity equation, with the density fluctuations that mostly seed the tempera- 
ture anisotropies. In a more detailed analysis of consistency, carried out in [69], a mild 
tension is reported between the QUaD data (particularly for Cf^) and the best-fitting 
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FIGURE 6. Cun-ent measurements of Cj^ (top) and Cf (bottom) from WMAP5 (magenta; [29]), 
QUaD (black; [63]), BOOMERanG (blue; [64, 65]), DASl (cyan; [66]), CAPMAP (green; [67]) and CBl 
(orange; [68]). The line is the best-fit flat, ACDM model to flie WMAP5 data. 



model to the WIVIAP data. The amplitude of the measured acoustic oscillations in C/ 
is rather higher for I > 500 than the best-fit model and this appears to be responsible for 
a high (probability to exceed of 7%) between the QUaD polarization data and the 
best-fit model. There will be significant improvements in measurements of the E-mode 
power spectra in the next three years, most notably from Planck. In addition, several of 
the new breed of polarimeters targeting 5-mode polarization have the angular resolution 
and sensitivity to measure E modes at high signal-to-noise to / ~ 2000, but their spectra 
will likely not be competitive with Planck due to their limited sky coverage. It will be 
interesting to see if the tension in the QUaD polarization data persists in these future 
datasets. 

4. COSMOLOGICAL PARAMETERS FROM THE CMB 

In this section we briefly discuss the dependencies of the CIVIB power spectra on the 
main cosmological parameters and the current measurements of the parameters via 
this route. Since the CIVIB observables are mostly a projection of conditions on the 
last-scattering surface, parameters can have an influence either through the acoustic 
physics of the pre-recombination plasma or through the angular diameter distance to 
last scattering, cIa - The latter controls the projection of linear distances at recombination 
to observed angular scales on the sky. We defer discussion of the effect of the primordial 
power spectra to Sec. 5.1. 
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The parameters that influence the perturbations up to the time of recombination are 
only the baryon and CDM densities - given that we know the CMB temperature rather 
precisely - in models with massless neutrinos. These densities set two length scales, 
the sound horizon rg and the damping scale and influence the relative heights of 
the acoustic peaks in the CMB power spectra (see below). The modulation of the peak 
heights largely determines the densities so that rg and K^^ become standard rulers at 
the time of last scattering. The angular scale of the anisotropics depends on the angles 
that these standard rulers subtend today, and hence on d-A- The main influence of the 
geometry of space, dark energy and light (but not massless) neutrinos on the observed 
anisotropies is through their effect on cIa- They also have a significant effect on the 
growth of structure at late times, but this only influences the linear CMB anisotropies 
through the late-time ISW effect on large scales where cosmic variance is large. 

To zero order, the primary CMB anisotropies and polarization determine the physical 
baryon and CDM densities, 0.\ih^ and Q.ch^, and cIa- Large-angle polarization further 
determines the optical depth to reionization. Models with the same values for these pa- 
rameters, and the same primordial power spectra, produce essentially the same CMB 
power spectra except on the largest scales [70]. This leads to a very precise geomet- 
ric degeneracy with CMB data alone that can only be broken with external data, e.g. 
the Hubble constant, luminosity distances from supemovae, and ^^(z) from the baryon 
acoustic oscillation (BAO) feature in galaxy clustering, or high-resolution observations 
of secondary effects in the CMB. A consequence of the geometric degeneracy is that 
the WMAP data alone is consistent with a model with a closed geometry and no dark 
energy, albeit with a very low Hubble constant (e.g. [71]). There are other approximate 
degeneracies in current CMB data, such as between the shape of the primordial pertur- 
bation spectrum (i.e. its spectral index, n^; see Sec. 5.1) and the matter density or the 
tensor-to-scalar ratio, r. 



4.1. IVIatter composition 

The dependence of the temperature and E-mode polarization power spectra on ^b^^ 
and Q.f^h^ is illustrated in Fig. 7. The main effects can be understood as follows. Consider 
increasing Q.\)h^ at fixed matter density Q.jnh^ = ^^b^^ + ^c^^- The increase in baryon 
inertia reduces the sound speed in the pre-recombination plasma by reducing its bulk 
modulus. As explained in Sec. 2.5, this shifts the mid-point of the acoustic oscillations 
to greater overdensities and enhances the height of the compressional peaks {n= 1,3,... 
for adiabatic initial conditions) in Cj . The sound horizon is also reduced, shifting the 
acoustic peaks in temperature and polarization to smaller scales (larger /). The increase 
in the number density of electrons in the plasma reduces the photon mean-free path, Ip, 
reducing the amount of diffusion damping and so increasing power on small scales. 
On larger scales, where tight-coupling still holds approximately, the reduction in Ip 
decreases the quadrupole anisotropy around recombination (the reduction in the plasma 
bulk velocity also contributes the same way) and the polarization is reduced. 

Consider instead increasing Q-mh^ at fixed Qb^^- At fixed dark energy density and 
curvature radius, the expansion rate would increase at all redshifts and dA would fall. 



Lecture notes on the physics of cosmic microwave background anisotropies 



March 30, 2009 



28 



500 1000 1500 2000 500 1000 1500 2000 
L I 

o 




500 1000 1500 2000 500 1000 1500 2000 
I I 



FIGURE 7. Variation of Cj (left), Cf (top right) and Cj^ (bottom right) as Hb/i^ is varied with fixed 
ilm/z^ (top) and as Q-mh^ is varied at fixed Hb/i^ (bottom). For the latter the angular diameter distance to 
last scattering has been held fixed. All models are flat. 



We have factored this effect out in Fig. 7 by reducing the dark energy density to 
keep d-A fixed. The main effects of increasing Q-mh^ are then to push matter-radiation 
equality back further in time and reduce the sound horizon and diffusion scale (since the 



Lecture notes on the physics of cosmic microwave background anisotropics 



March 30, 2009 



29 



conformal age of the universe is reduced). Since the residual radiation is then lessened at 
recombination, the early-ISW effect and thus the first peak in Cf are also reduced. The 
resonant driving effect that enhances the amplitude of the acoustic oscillations for modes 
that enter the sound horizon during radiation domination is also limited to higher peak 
order. Together, these effects lead to the relative enhancement of the third and higher 
peaks over the first two that can be seen in Cf in Fig. 7. The reduction in the sound 
horizon shifts the temperature and polarization peaks to larger /. Decreasing diffusion 
reduces the polarization on large scales where tight-coupling holds, as does the lessened 
impact of resonant driving. 

These effects on the morphology of the acoustic peaks have allowed accurate mea- 
surements of the matter and baryon densities from the CMB. Five years of WMAP data 
alone give t^b/z^ = 0.02273 ±0.00062 (i.e. 3% uncertainty) and Qc^?^ = 0.1099 ±0.0062 
(and so Qmh^ = 0.1326 ± 0.0063; 6% uncertainty) in flat ACDM models [9]. These 
numbers should improve to sub-percent levels with better measurements of the third and 
higher peaks with the Planck data in the future [4]. 

Neutrino oscillations imply that at least two flavours are massive (see [72] for a review 
in the cosmological context). The minimum mass of the heaviest eigenstate is ?a 0.05 eV. 
At this minimum level, neutrinos are still relativistic at recombination and their effect on 
the dynamics of the primordial plasma are indistinguishable from massless neutrinos. A 
neutrino is non-relativistic at recombination if its mass exceeds ^ 0.6 eV. At this mass, 
the eigenstates must be very nearly degenerate and the total (summed) neutrino mass is 
therefore Ly'^v ~ 1.8eV. This sets a rough limit on the determination of masses from 
the primary CMB anisotropics and is close to the limit Y.v'^v < l-3eV from the five- 
year WMAP data alone in flat, A models [44]. Even with minimal masses, two of the 
neutrino eigenstates are non-relativistic today and therefore affect and suppress the 
late-time growth of structure on scales small enough that the neutrinos cannot cluster. 
The former effect can be used to improve constraints on Y,v f^v from the CMB if external 
data is used to break the geometric degeneracy. In WMAP5, the constraint on Y.v'^v 
improves to 0.67eV in flat, A models by including supernovae and BAO distance data. 
The effect of sub-eV neutrinos on structure formation affects the lensing of the CMB 
by large-scale structure, making CMB lensing a promising technique for determining 
masses (see Sec. 6.1.2). 



4.2. Geometry and dark energy 

Generalising our discussion of flat models in Sec. 2.5 to non-flat models, the positions 
of the acoustic peaks in models with adiabatic initial conditions are asymptotically at 
/ = mtdA/rs{r\*) since curvature is never dynamically important before recombination. 
In flat models, note that d-A reduces to x*^ the comoving distance to last scattering. With 
the sound horizon determined from the peak morphology, the CMB alone can be used 
to measure dA'. constraints from WMAP5 give ^4 = 14. 1 ± 0.2 Gpc in their most general 
open models with dynamical dark energy [44]. 

The dependence of Cj on the curvature fraction, Q.k (i.e. the ratio of the current 
Hubble radius to curvature radius and defined to be positive/negative for open/closed 
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models), and dark energy density in ACDM models is illustrated in Fig. 8. Both param- 
eters principally affect the anisotropics through and so simply shift the peaks. The 
remaining effect on large scales is due to the late-time ISW effect and, in closed models, 
from mode quantisation. Disentangling the effects of dark energy and curvature from 
the linear CMB anisotropics requires external data to break the geometric degeneracy. 
Combining WMAP5 with BAO and supernova data gives -0.0178 < ^^a: < 0.0066 in A 
models, fully consistent with flatness. (Only one external dataset is required in this case; 
BAO is the most constraining.) The geometric degeneracy also limits constraints from 
the CMB alone on more complex models with dynamical dark energy, even for a flat uni- 
verse; a detailed discussion is given in [44]. The current constraints combining WMAP 
with both BAO and supernova data are consistent with flatness and non-dynamical dark 
energy: the equation of state parameter w is consistent with —1 at the 15% level. 



5. CONSTRAINING EARLY-UNIVERSE PHYSICS WITH THE 

CMB 

5.1. Inflation and the origin of structure 

So far, we have not discussed the origin of the primordial perturbation which pro- 
vided the seeds for cosmological structure formation under the action of gravitational 
instability. The study of this question has the potential to expose deep connections be- 
tween cosmology and physics at immensely high energies which are forever beyond the 
reach of earth-bound particle accelerators. Inflation [73, 74, 75, 76, 77, 78], an epoch 
in which the expansion of the universe is accelerating, solves a number of puzzles asso- 
ciated with the standard big bang cosmology. During a phase of accelerated expansion, 
(the physical Hubble radius) remains almost fixed but the physical separation of 
particles initially in causal contact grows exponentially. The result is that regions today 
separated by cosmological distances were actually in causal contact before/during infla- 
tion. At that time, these regions were given the necessary initial conditions, smoothness 
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and small perturbations about smoothness, we observe today. 

The big bang puzzles (e.g. horizon problem, flatness problem, homogeneity problem 
and monopole problem) are not strict inconsistencies in the standard big bang theory, but 
must instead be assumed as extreme fine-tuning of the initial conditions, thus severely 
restricting the predictive power of the theory. Inflation provides a dynamical mechanism 
to solve the big bang problems. At the heart of the inflationary paradigm is an elegantly 
simple idea: invert the behaviour of the comoving Hubble radius i.e. make it decrease 
sufficiently in the very early universe. It is easy to see from the second Friedmann equa- 
tion that the three equivalent conditions necessary for inflation are decreasing comoving 
Hubble radius, accelerated expansion, and violation of the strong energy condition (i.e. 
negative pressure): 

The requirement for negative pressure means that inflationary dynamics cannot be 
achieved via normal matter or radiation! Inflation can be phenomenologically described 
by a scalar field with special dynamics, evolving in a self-interaction potential. Although 

no fundamental scalar field has yet been detected in experiments, there are fortunately 
plenty of such fields in theories beyond the standard model of particle physics. Can a 
scalar field have p -I- 3p < 0? The equation of state for a homogeneous scalar field 4>, 

shows that a scalar field can lead to accelerated expansion (w < —1/3) if the potential 
energy V dominates over the kinetic energy {dt^)^/2. The equation of motion of such a 
homogeneous scalar field is 

d^^ + 3Hdt^ + V^ = (75) 

where V<d = dy/d4>. 

In the following discussion, for notational convenience we will set the reduced Planck 
mass to unity, i.e. SkG = Mp^^ = 1. As we have seen, inflation occurs if the field is 
evolving slow enough that the potential energy dominates over the kinetic energy, and in 
order to solve the classical big bang problems, the second time derivative of 4> is small 
enough to allow this slow-roll condition to be maintained for a sufficient period. Thus, 
successful inflation requires 

{dt^f < y(4>) (76) 
\d^^\ < \3Hdt^\,\V<p\, (77) 

where the latter condition means the motion is friction dominated, i.e. 3Hdf(^ ~ — 
from Eq. (75). Satisfying these conditions requires the smallness of two dimensionless 
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quantities known as potential slow-roll parameters 




2 



(78) 



V 



(79) 



In the slow -roll regime, Cy, |t]v| <^ 1, with background evolution 



d4> 



y(4>) K, const., 
37f ' 



(80) 



(81) 



and the spacetime is approximately de Sitter, with the scale factor evolving as a{t) ~ e 
where the Hubble parameter H is approximately constant. 

Inflation ends when the slow-roll conditions are violated: ey(^end) ~ 1- The number 
of e-folds before inflation ends is 



Roughly 60 e-folds of inflation must occur for the classic big bang puzzles to be 
resolved. The universe must then reheat, the energy density in the scalar field being 
converted to radiation to start off the radiation-dominated era. 

Besides solving the big bang puzzles, the decreasing comoving horizon during in- 
flation is the key feature required for the quantum generation of cosmological pertur- 
bations. During inflation, quantum fluctuations are generated on sub-Hubble scales and 
are then stretched out of the Hubble radius by the accelerated expansion. In other words, 
the superluminal expansion stretches the perturbations to apparently acausal distances. 
They become classical superhorizon density perturbations which re-enter the Hubble 
radius in the subsequent non-accelerating evolution and then undergo gravitational col- 
lapse to form the large-scale structure in the universe. 

To analyse fluctuations during inflation, the inflaton field is split into a homogeneous 
background 4>(r) and a spatially varying perturbation 54>(r,x): 



Perturbations of the inflaton field value satisfy the equation of motion of a harmonic 
oscillator with time-dependent mass. The quantum treatment of inflaton perturbations 
therefore parallels the quantum treatment of a collection of one-dimensional harmonic 
oscillators. Just as zero-point fluctuations of a harmonic oscillator induce a non-zero 
variance for the oscillation amplitude {x^), the quantum fluctuations of a light scalar 
field^ during inflation induce a non-zero variance for the inflaton perturbations [79]. The 



^ Such a field has (effective) mass y<Do <C which is equivalent to |tJv| <C 1. 
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(82) 



^>(?,x) = ^)(?) + 54>(?,x). 



(83) 



variance of these fluctuations in Fourier space, i.e. the power spectrum, is 



(84) 



where the right-hand side is evaluated when the mode exits the Hubble radius. Our 

Fourier conventions follow those in Sec. 2.4.1. Technically, the field fluctuation here 
is defined on hypersurfaces with zero intrinsic curvature and, in this gauge, the result 
is the same as if metric perturbations (i.e. the back-reaction of 54> on the spacetime 
geometry) were ignored. 

A pedagogical review of the generation of density (scalar) perturbations is beyond the 
scope of this article and we direct the interested reader to the comprehensive treatment 
in Ref. [80]. The quantity we wish to compute is the comoving curvature perturbation 
since this is conserved from Hubble exit during inflation to Hubble re-entry during 
the standard radiation or matter- dominated epochs. With ^ we can reliably compute the 
primordial fluctuation (in single-field models of inflation) in the radiation era without 
needing to model the dynamics of the reheating process. During inflation, comoving 
hypersurfaces have the property that they coincide with the hypersurfaces over which 
the (total) inflaton 4> is homogeneous. These hypersurfaces are not the same as the 
zero-curvature surfaces on which Eq. (84) holds - there is a time delay between them, 
5t — — 54>/ (9f<l>, such that the evolution of the background <l> in this time compensates 
for the perturbation 54> to give a smooth total The differential background expansion 
during this time delay means that the intrinsic curvature of the comoving hypersurfaces 
is simply M = —Hd^/dt^. In this simple class of inflation models, the power spectrum 
of ^ is thus [81] 



where the second equality uses the slow-roll approximations (80) and (81). The effi- 
ciency of this mechanism for producing cosmological curvature perturbations depends 
on both the height of the potential, which determines the expansion rate and hence size 
of and its slope, which enters through the conversion of inflaton fluctuations to time 
delays and so curvature. 

For gravitational waves, the two polarization modes of the metric perturbations hij 
(see Eq. 37) satisfy the same equation of motion as the perturbations of a massless scalar 
field. The power spectrum of hij is therefore [77, 78, 81] 



where the second equality uses the slow-roll approximation. Our Fourier conventions 
for hij are such that the real-space variance of the metric fluctuation is 



(85) 




(86) 




(87) 
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Note that the tensor power spectrum depends only on the expansion rate during inflation. 
The amplitude of the tensor power spectrum relative to the curvature spectrum defines 
the tensor-to-scalar ratio 



= 16ey. (88) 



Here, k^, is some (arbitrary) reference scale. As we will explain below, the value of r 
is of fundamental importance in the quest to understand the microscopic origin of the 
inflationary dynamics. 

The power spectra in Eqs. (85) and (86) are to be evaluated when a fluctuation 
with (physical) wavenumber k/a exits the Hubble radius H^. Different scales exit at 
different times when the inflationary potential has slightly different values. This 
leads to a small (but computable) scale-dependence of the primordial power spectra. In 
terms of an empirical parameterisation of this weak scale-dependence, the scalar power 
spectrum can be written in terms of an amplitude As and a spectral index as 

j^^{k)=Mk/k^r^-\ (89) 

both measured at the reference scale k^,. Similarly, its tensor counterpart can be written 
as 

^h{k)=Mk/k^y'. (90) 

The deviation from scale invariance (n^ = 1, = 0) of these power spectra is sensitive 
to the shape of the inflaton potential: 

„,-l = 2%i-3fi?f , (91) 



and 



V 



n, = -4[^]\ (92) 



The primordial perturbation spectra are rather directly imprinted in the observable CMB 
angular power spectra. For example, for scalar perturbations, the dominant Sachs- Wolfe 
contribution to the anisotropics gives Q °^ ^/^(//J^). where the /-dependent proportion- 
ality depends on the processing of the primordial perturbations by gravity and acoustic 
physics. Measurements of the scale-dependence of the primordial power spectra there- 
fore have the power to probe the shape of the inflaton potential V(^»). Note the consis- 
tency relation, r = — 8nf, in the simple single-field models we are considering so far. 

Two important pieces of information about inflation would follow from a detection of 
A;: (1) the energy scale of inflation; and (2) the field variation of the inflaton. The energy 

scale of inflation V^^^ is proportional to A^^. Using the measured value for A^, this can 
be expressed in terms of the tensor-to-scalar ratio on CMB scales, r, as 

r \i/4 



y^/'^ = 1.06X 10^^GeV(^^) . (93) 
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To date, high-energy physicists only have two indirect clues about physics at this scale: 
the apparent unification of gauge couplings, and experimental lower bounds on the 
proton lifetime, and these energy scales are forever beyond the reach of earth-bound 
particle accelerators. Thus, if a primordial tensor mode were to be detected, we would 
be presented with a unique opportunity to use the universe as a high-energy physics 
laboratory. 

From manipulations of Eq. (88) one may derive [82] the following relation between 
the tensor-to-scalar ratio r and the distance in inflaton field space between the end of 
inflation and the point when the scales of CMB fluctuations were created: 



Here, we have reinstated the reduced Planck mass. A large tensor amplitude, r > 
0.01, therefore correlates with a super-Planckian field variation during inflation. Super- 
Planckian field excursions have interesting theoretical implications [83]: to control the 
shape of the inflaton potential over a super-Planckian range requires the existence of an 
approximate shift symmetry in the ultraviolet (UV) limit of the underlying particle the- 
ory for the inflaton, ^> ^ + const. In string theory it has only recently become possible 
to construct controlled large-field inflation models with approximate shift symmetries in 
the UV [84, 85, 86]. 



5.2. Observable predictions and current observational constraints 

The inflationary proposal requires a huge extrapolation of the known laws of physics. 
In the absence of a complete theory, a phenomenological approach has been commonly 
employed, where an effective potential is postulated. Ultimately, V{^) has to be 
derived from a fundamental theory, and significant progress in implementing inflation 
in string theory has been made in recent years [86]. However, while it is challenging to 
understand the origin of inflation from a particle physics point of view, it is also a great 
opportunity to learn about ultra-high-energy physics from cosmological observations. 

The simplest inflationary scenarios consist of a single light scalar field with a canon- 
ical kinetic term, (V4>)^/2, in its action. They predict the following observable charac- 
teristics. 

1. Flat geometry, i.e. the observable universe should have no spatial curvature. As 

we have seen, flatness has been verified at the 1% level by the location, or, better, 
separation, of the CMB acoustic peaks combined with some low-redshift distance 
information. 

2. Gaussianity, i.e. the primordial perturbations should correspond to Gaussian ran- 
dom variables to a very high precision. 

3. Scale-invariance, i.e. to a first approximation, there should be equal power at all 
length-scales in the perturbation spectrum, without being skewed towards high or 
low wavenumbers. In terms of the parameterisation (89) and (90) this corresponds 
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to = 1 and — 0. However, small deviations from scale-invariance are also a 
typical signature of inflationary models and tell us about the dynamics of inflation. 

4. Adiabaticity, i.e. after reheating, there are no perturbations in the relative number 
densities of different species on super-Hubble scales (so no isocurvature modes). 
This follows from the assumption that only a single field is important during 
inflation. Constraints on isocurvature modes were discussed in Sec. 2.7.4 where 
we noted that current data shows no evidence for a non-adiabatic component of the 
primordial perturbation. 

5. Super-Hubble fluctuations, i.e. there exist correlations between anisotropics on 
scales larger than the apparent causal horizon, beyond which two points could not 
have exchanged information at light-speed during the history of a non-inflationary 
universe. This corresponds to angular separations on the sky larger than ~ 2°. 

6. Primordial gravitational waves, which give rise to temperature and polarization 
anisotropics as described above. These tensor modes must exist; however, their 
predicted amplitude can vary by many orders of magnitude depending on the 
underlying microphysical mechanism implementing inflation. 

It is beyond the scope of these lecture notes to describe in algorithmic detail how 

CMB observations can be used to constrain inflationary physics; the interested reader is 
referred to Ref. [87] for a beginner-level introduction and a guide to the primary liter- 
ature. The tightest constraints on inflationary physics currently come from the WMAP 
data in combination with complementary cosmological data. We now summarise some 
of the main consequences of the five-year data for inflation where these have not been 
covered elsewhere in these notes. 



5.2.1. Super-horizon correlations 

The WMAP detection of an (anti-)correlation between CMB temperature and po- 
larization fluctuations at angular separations 5° > > 1° (corresponding to the TE 
anti-correlation seen for multipoles / ~ 50-150 in Fig. 6) is a distinctive signature of 
adiabatic fluctuations on super-horizon scales at the epoch of recombination, confirming 
a fundamental prediction of the inflationary paradigm [57, 88, 89]. Inflation, in which 
microscopic quantum fluctuations are stretched to super-Hubble scales, is the most com- 
pelling causal mechanism for generating such apparently-acausal perturbations. The ob- 
served TE anti-correlation on large scales is remarkable qualitative evidence for this 
basic mechanism. 



5.2.2. Measurements of the scalar spectrum 

Komatsu et al. [44] recently used the WMAP five-year temperature and polarization 
data, combined with the luminosity distance data of Type-la supemovae (SN) at redshifts 
z < 1.7 [90] and BAO data at redshifts z = 0.2 and 0.35 [91], to put constraints on the 
shape of the primordial power spectra (see Table 1). 
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TABLE 1. Five-year WMAP constraints on 
the primordial power spectra in the power-law 
parameterisation [44]. 





5-year WMAP 


WMAP+BAO+SN 


n.s- 


963+0 014 


0.960±0.013 


n,v 


0.986 ±0.022 


0.970 ±0.01 5 


/• 


<0.43 


< 0.22 


ris 

CCs 


1 031+0.054 

-0.037 ±0.028 


^•^^ ' -0.043 

-0.028±0.020 


ris 
r 


1 087+0 072 

<0.58 


1 089+00™ 
<0.55 


CCs 


-0.050 ±0.034 


-0.053 ±0.028 



The WMAP analysis employed the standard power-law parameterisation of the 
scalar power spectrum given in Eq. (89). The amplitude of scalar fluctuations at 
h = 0.002 Mpc-i is measured to be As{h) = (2.445 ±0.096) x 10"'^. From both nu- 
merical simulations and analytical estimates one finds that with this initial amplitude of 
density fluctuations, there is sufficient time for gravity to form the large-scale structures 
that we observe today. Assuming no tensor perturbations (r = 0) the scale-dependence 
of the power spectrum is 

n, = 0.960 ±0.013 (r = 0). (95) 

The scale-invariant Harrison-Zel'dovich-Peebles spectrum, = 1, is 3.1a away from 
the mean of the likelihood. Including the possibility of a non-zero r in the parameter 
estimation, the marginalised constraint on Ug becomes 

n, = 0.970 ±0.015 (r^O). (96) 

The (slight) worsening of the constraint simply means that there is a degeneracy between 
Us and r, and the current data cannot simultaneously measure them independently. Since 
adding tensor perturbations enhances the large-scale power, it can be offset by increasing 
rig. 



5.2.3. Constraints on the spectrum of tensor perturbations 

At current sensitivities, constraints on the amplitude of tensor perturbations are driven 
by the temperature data. In the future, this situation should change with weight shifting 
to the 5-mode of polarization, as discussed later in Sec. 5.3.2. With five years of WMAP 
data alone, one finds the following upper limit on the tensor-to-scalar ratio: 

r<0.43 (95%C.L.). (97) 

This improves to r < 0.22 (95% C.L.) if SN and BAO data are included. The latter 
datasets add only geometric information and so are not influenced by tensor perturba- 



Lecture notes on the physics of cosmic microwave background anisotropies 



March 30, 2009 



38 



tions. However, what they do bring is a better determination of the matter density which, 
through degeneracies, leads to an improvement in and hence r. The Umit r < 0.22 is 
close to the fundamental limit, set by cosmic variance, that can be achieved with temper- 
ature data alone. Upcoming polarization experiments with sensitivity to 5-modes should 
improve on this limit by an order of magnitude or more. 



5.2.4. Constraints on inflationary models 

At the present stage we are still testing basic aspects of the inflationary mechanism 
rather than details of its specific implementation. However, constraints on specific in- 
flationary models have been obtained (see Fig. 9): inflationary models predicting a blue 
spectrum (n^ > 1) are now virtually ruled out; this includes models of hybrid inflation 
like V{(^) = Vq[1 -|-m^4>^]. Assuming that the tensor amplitude is small, the measure- 
ment Us < I implies a constraint on the curvature of the inflaton potential, V" < 0. Fi- 
nally, models predicting a very large tensor amplitude are ruled out, e.g. V{^) = X^"^. In 
summary, many popular models are still allowed by the data, but an increasing number 
of models are on the verge of being tested seriously. 



5.3. Future prospects for constraining early-universe physics with the 

CMB 

There are exciting prospects for improving our understanding of early-universe 
physics with new CMB data that will be available in the next decade. The Planck satel- 
lite [4] will measure the temperature power spectrum over a large range of scales with 
unprecedented accuracy and resolution. In addition, it will provide greatly improved 
constraints on E- and 5-mode polarization. The Planck data will be supplemented by 
many sub-orbital experiments with a special focus on measurements of the small-scale 
temperature fluctuations and/or the polarization power spectra. Finally, plans are being 
made for a next-generation satellite dedicated to the measurement of CMB polariza- 
tion [83, 99]. Such an experiment proposes to improve the sensitivity to 5-modes by 
almost two orders of magnitude over the current situation. As we now describe, the 
combination of this wealth of data will allow detailed tests of the physics of the early 
universe. 



5.3.1. Scalar perturbation spectrum on small scales 

The shape of the spectrum of primordial density fluctuations can be an important 
diagnostic of the inflationary dynamics. Equation (89) is, of course, only a simple 
parameterisation of the power spectrum in terms of an amplitude and a spectral index 
Hs, both defined at the pivot scale A;^. This power-law parameterisation may be refined by 
allowing a non-zero running of the spectral index, a, = dus/dlnk = d^hi^^/d{\nk)^. 
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FIGURE 9. Constraints on single-field slow -roll models in the n^-r plane [83]. The value of r deter- 
mines whether the model involves large or small inflaton field variations. The sign of Hs — 1 classifies the 
scalar spectrum as red or blue. The combination 2tjv = Hj — 1 + 3r/8 determines whether the curvature of 
the potential was positive (rjy > 0) or negative (r\Y < 0) when the observable universe exited the Hubble 
radius. Shown also are the WMAP five-year constraints on and r [44], as well as the predictions of a 
few representative models of single-field slow-roll inflation: chaotic inflation - Xp 4>'', for general p (thin 
solid line) and for p — 4,3,2, 1,2/3 (filled circles) - models of which with p = \ [92] and p — 2/3 [84] 
have recently been obtained in string theory; natural inflation - V{,)[1 — cos(<I>//x)] (solid line); hill-top 
inflation - Vo[l — (<!>/ /i)^] + . . . (solid line); very small-field inflation - models of inflation with a very 
small tensor amplitude, r <^ 10^^ (green bar), e.g. brane inflation [93, 94, 95, 96], Kahler inflation [97], 
and racetrack inflation [98], which often arise in string theory. 



i.e. by defining 

/ J \ )ij(^,t)-l+aj(^v.)ln(fc//:*)/2 

^.^{k)=A,{K){-^\ . (98) 

Generically, slow -roll inflation predicts that the running should be a small effect, second 
order in the slow-roll parameters, and so la,! ~ 0(0.001). At present, the data is insuf- 
ficiently precise to detect Us at that level. The current constraints are rather dependent 
on the combination of CIVIB data used, and the details of the treatment of contami- 
nants at high multipoles such as emission from extra-Galactic sources and the Sunyaev- 
Zel'dovich effect in galaxy clusters. Full details can be found in [30, 44]; the basic 
result is that no running (a, = 0) is consistent with various data combinations at the 
95% confidence level but the mean values of Us are persistently negative. The key to 
improving constraints on running is to increase the lever arm in k. Future small-scale 
CIVIB measurements, possibly combined with large-scale structure surveys, should im- 
prove constraints considerably though probably not to the level where the prediction of 
minimal inflation becomes detectable. Forecasts for Planck suggest la (marginalised) 
errors of 0.005 for a, [4]. 
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In the case of slow-roll inflation, a definitive measurement of runnning of the spectral 
index, Ug 7^ 0, is a signal that ^h, the third Hubble slow-roll parameter (defined in 
analogy to the first two potential slow-roll parameters discussed previously),^ 



(99) 



played a significant role in the dynamics of the inflaton [100] as the CMB scales exited 
the Hubble radius. In the context of slow-roll inflation it would be hard to explain a 
value of Ug that is either negative or positive but much larger in magnitude than 0.001. 
The consequences for the physics of inflation differ depending on whether the running 
is negative or positive, and both options would dramatically complicate the theoretical 
understanding of inflation. 

A large negative running implies that was (relatively) large and positive as the 
cosmological perturbations were laid down (a^ ^ —2^h in this limit). It can be shown 
that > generally hastens the end of inflation (relative to — 0), provided the 
higher order slow-roll parameters can be ignored. With these assumptions, we find a 
tight constraint on if we are to avoid a premature end to slow roll, with inflation 
terminating soon after the observable scales exit the Hubble radius [101, 102, 103, 
104]. Thus, a definitive observation of a large negative running would imply that any 
inflationary phase requires higher order slow-roll parameters to become important after 
the observable scales leave the Hubble radius [101, 102, 105, 106], multiple fields 
yielding complicated spectra [107], two or more bursts of inflation [108, 109, 110, 
111, 112], or a temporary breakdown of slow roll (such as a feature in the potential) 
[89, 113, 114, 115, 116, 117, 118]. There have also been proposals for explicit models 
where a large, transient negative running occurs as a "stringy" signature in multi-throat 
brane inflation [119, 120, 121]. 

The current cosmological data disfavour inflationary models with a blue tilt, > 1 
[9, 44]; however, a significant parameter space is still allowed where ^5 < 1 but with a 
large positive running (implying a large negative ^h), which would lead to a strongly 
blue-tilted spectrum on scales below those of cosmological interest [122]. Again under 
the hypothesis that this parameterisation can be extrapolated to the end of inflation, we 
find a class of solutions where e — > as remains finite, and the field rolls towards a 
minimum with a substantial vacuum energy. The perturbation spectrum grows at small 
scales, possibly diverges, and can lead to over-production of primordial black holes 
[122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133], or even the onset of eternal 
inflation [122, 134, 133]. Consequently, if we measure a large positive running we will 
again conclude that any inflationary phase is not described within the single field, slow- 
roll formalism, or that higher order terms in slow roll are important. 



The primes here denote derivatives with respect to 3>. 



Lecture notes on the physics of cosmic microwave background anisotropies 



March 30, 2009 



41 




10 100 1000 



I 

FIGURE 10. Cun-ent 95% upper limits on Cf from WMAP5 (magenta; [29]), QUaD (black; [63]), 
BOOMERanG (blue; [64, 65]), DASI (cyan; [66]), CAPMAP (green; [67]), CBI (orange; [68]) and 
POLAR (dashed magenta; [135]). The solid line is the contribution from gravitational waves for a tensor- 
to-scalar ratio r = 0.2. The dashed line is the B-mode power induced by weak gravitational lensing. 



5.3.2. Tensor perturbations 

The WJMAP upper limit on the amplitude of primordial tensor fluctuations, r < 0.22, 
is only beginning to be a serious constraint on inflationary model-building. As noted 
above, current constraints are driven by the large-angle temperature data, but this route 
is fundamentally limited by the cosmic variance of the dominant scalar perturbations. 
The same is true for ii-mode polarization but not for 5-mode polarization since this is 
not generated by scalar perturbations in linear theory. Improved constraints on 5-modes 
of CIVIB polarization will therefore drive future limits on the tensor amplitude. The 
quest for a 5-mode detection is one of the most exciting developments in observational 
cosmology. 

Current limits on 5-mode power are compared in Fig. 10 to the power expected 
for r = 0.2. (The secondary contribution from conversion of E modes to B by weak 
gravitational lensing is also shown; see Sec. 6.1.2.) The primordial signal is very weak, 
with r.m.s. = 0.lS^/rj0?2iJ.K, but there are no fundamental technical or cosmological 
obstacles to achieving r ~ 0.001, although it calls for arrays of thousands of detectors. 
However, it is still unclear whether polarized Galactic emission, particularly thermal dust 
emission, is sufficiently well behaved to allow its removal down to this level on the basis 
of multi-frequency data. There are two characteristic scales for attempting detection. 
The reionization signal at / < 10, and the signal from recombination at Z < 100. Where 
most information lies depends on the value of r, the instrument sensitivity, the scale- 
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dependence of Galactic foregrounds, and the size of the survey (i.e. whether the large- 
scale reionization modes are accessible). Realistic forecasts for Planck suggest r = 0.05 
may be achievable, and the constraint is driven by the reionization signal [136]. In 
contrast, most of the next generation of ground-based and balloon experiments will 
survey smaller regions of the sky in areas of (known) low foreground contamination, 
and their target is the signal from recombination [137, 138, 139, 140, 141, 142]. Tensor- 
to-scalar ratios r ~ 0.01 may be achievable with these surveys. Looking further ahead, 
r ~ lO""' with a full-sky survey from a dedicated polarization satellite seems achievable 
with current technology. Measuring tensors at that level of sensitivity would mark a 
qualitative shift in our ability to test the inflationary paradigm. If tensors are not seen at 
that level, all large-field models of inflation are ruled out. On the other hand, a 5-mode 
detection would be an extraordinary discovery and a "smoking gun" for inflation. 

A natural target sensitivity for most future fi-mode experiments is ~ 5 /xK-arcmin in 
the polarization maps. At this level, the "noise" due to the B modes produced by gravi- 
tational lensing is comparable to the instrument noise [54, 143]. For reference, this is 10 
times the sensitivity of the nominal Planck survey. Pushing for further improvements in 
sensitivity is then only worthwhile if it improves rejection of foreground contamination, 
or if the lensing contamination can be removed coherently (i.e. at the level of the maps), 
for example with quadratic reconstruction techniques [144]. Without reconstruction, the 
(3a) upper limits on r set by lensing, assuming r = 0, are 1 x 10^^ from the reionization 
signal and 8 x 10^ from recombination alone (i.e. / > 10). In practice, astrophysical 
foregrounds and instrumental systematic effects will surely be more of a limitation than 
lensing noise. Foregrounds will be particularly troubling for the signal from reioniza- 
tion, with estimates suggesting they may limit detections to r ~ 0.01 (99% C.L.) from 
/ < 10 data alone given the frequency coverage and sensitivity being considered for a 
future satellite mission [145, 146]. 



5.3.3. Non-Gaussianity 

So far, we have assumed that the primordial fluctuations have a Gaussian distribution. 
Indeed, this is a fundamental prediction of slow-roll inflation and the current observa- 
tional limits confirm this at the 0.1% level.^ However, as we will now discuss, a small 
degree of primordial non-Gaussianity can be a crucial probe of the inflationary dynam- 
ics. 

Non-Gaussianity is a sensitive probe of inflationary physics that is difficult to access 
by any other means. Specifically, it is a measure of inflaton interactions. To allow for 
slow-roll inflation the inflaton field is necessarily weakly interacting [y(4>) is very 
flat] and the non-Gaussianity is predicted to be small [147]. However, going beyond 
the single-field slow -roll paradigm, non-trivial kinetic terms (derivative self-interactions 
in the inflaton action), the presence of more than one light field during inflation, the 
(temporary) violation of slow roll, and a non-adiabatic initial vacuum state for the 
inflaton may lead to large, observationally distinct non-Gaussianity (e.g. [148, 149]). 



' By this measure, non-Gaussianity has been constrained more accurately than curvature. 
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Gaussian fluctuations are characterised completely by their two-point correlation 
function, or, equivalently, their power spectrum, ^ag{k). Non-Gaussianity is therefore 
measured by considering the connected higher-order correlation functions, i.e. the part 
that remains after all possible contractions are subtracted. The connected n-point func- 
tions with n> 2 vanish for a Gaussian random field and so the leading non-Gaussian 
effect is usually given by the three-point function. To calculate the non-Gaussianity from 
inflation, generally one must go beyond linear perturbation theory, expanding the action 
to third order to capture all cubic interactions between perturbed quantities [147]. The 
primordial 3-point function from inflation is conveniently calculated for quantities C, 
and Jij which are the non-linear generalisations of the linear-theory comoving-gauge 
curvature perturbation, and hij, the linear- theory gravitational wave amplitude, respec- 
tively.^ Both C, and jij are conserved on super-Hubble scales in the simplest inflation 
models. The bispectrum for C, is defined by 

(C(ki)C(k2)C(k3))=5^(/:i,/:2,fc3)5(ki+k2 + k3). (101) 

Three other bispectra with different combinations of C, and jij can also be constructed. 
Note that the bispectrum depends only on the magnitudes of the three wavevectors as a 
consequence of translational, rotational and parity invariance. 

The bispectrum is measured by sampling triangles in Fourier space. Much physical 
information is contained in the momentum dependence or the shape of the bispectrum. 
Ref. [150] presents visualisations of the full scale/shape dependence of the bispectra 
for a variety of inflationary scenarios in which significant (observable) amounts of non- 
Gaussianity could be potentially produced. While these have complex shapes, fortu- 
itously, different physical mechanisms that produce significant non-Gaussianity result 
in bispectrum signals that peak on different triangular configurations (e.g. [48, 149, 151, 
152, 153, 154]): 

• multi-field models peak on squeezed triangles (^3 <S ^2 ~ ^1); 

• models with non-canonical kinetic terms peak on equilateral triangles {k\ ~ /:2 ~ 

• models with non-adiabatic initial vacua peak on flattened/folded triangles (% ~ 
k2 ~ 2^1); and 

• non-slow-roll models peak on more complex (model dependent) configurations 
requiring a matched template for their analysis. 

Thus, configuration-dependent studies of non-Gaussianity may become a powerful 
probe of ultra-high-energy physics and inflation. 

Given a primordial bispectrum, how do we compute the observable bispectrum, ^/j/^/j, 
in the CMB (see Eq. 8)? In theories with weak primordial non-Gaussianity, such as 



° These quantities are defined so that the induced line element on comoving hypersurfaces, over which 
53> = 0, is [147] 

^g,;dx W = a^e^^e'^'' . (1 00) 

Here, jij is trace- free and transverse, difj = 0, so that dete'^^ = 1 . Note that the true curvature perturbation 
of comoving hypersurfaces differs from ^ by terms quadratic in C, and jij. 
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single-field inflation, this is a difficult task since the non-linear perturbation evolution 
and radiative transfer will induce additional non-Gaussianity in the CMB that exceeds 
(and distorts) the primordial contribution. Signatures of specific non-linear effects have 
been computed, but a full calculation is still lacking (see [155, 156, 157, 158, 159, 
160, 161] for some recent progress towards this goal). Things are more straightforward, 
though, if the initial level of non-Gaussianity is large enough, since then one can con- 
tinue to employ linear perturbation theory to compute the higher-order statistical prop- 
erties of the CMB, replacing the bispectrum of the linear-theory with 5^ (^i , ^2, ^3)- 

Generalising Eq. (22) to include all sources of anisotropy, the multipoles of the 
temperature anisotropy, T/^, can be written in the form 

^^gK^)^(k)C(k)^*"^ (102) 

in linear theory with adiabatic perturbations. Here, gi{k) is the temperature transfer 
function which linearly relates the observed CMB anisotropy to the primordial curvature 
perturbation. The observed CMB bispectrum then becomes 



_ ; (2/i + l)(2/2+l)(2/3 + l) //i /2 /3 



where the reduced bispectrum &/J/2/3 is [162] 

bhhh = I ^'dr J n [kfdkigi,ik)ji,ikir)] B^ihMM) ■ (104) 



As a simple and well-studied example, we consider the local model of primordial 
non-Gaussianity in which ^ is a sum of a Gaussian piece and the square of a Gaussian 
in real-space: 

C(x) = Cg(x) + ^/nl [caw - (C^)] • (105) 

The factor of 3/5 is conventional and is chosen so that — /nl parameterises the quadratic 
contribution to the gravitational potential ^ in matter domination (for I/nlI ^ 1 so fur- 
ther non-linear evolution can be ignored).^ Single-field, slow-roll inflation with canoni- 
cal kinetic term approximates to the local model, with a very small primordial /nl of the 
order of the slow-roll parameters, i.e. ~ 0(0.01) [147] (non-linear effects - e.g. in the 
relation between the gravitational potential and C - increase the effective observable /nl 
to 0(1)). Large local non-Gaussianity can be generated in models with multiple fields, 
such as the curvaton model. The bispectrum in the local model evaluates to 

Bi^ihMM) = (^^^M^^M+P^rms^ . (106) 



^ Our sign choices and normalisation are consistent with [44]. Note that on large scales, the temperature 
anisotropics are skewed negative by a positive /nl- 
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Significantly, this is factorisable which makes the integration in Eq. (104) easily 
tractable. Numerical examples of the reduced bispectrum in the local model can be 
found in Ref. [163]; it displays similar acoustic oscillations to the power spectrum. The 
best observational constraints on the local form of non-Gaussianity are from the 3-point 
function of the WMAP5 data. With an optimal weighting of the data, Ref. [164] find 
no evidence for non-zero /nl: —4 < /nl < 80 (95% C.L.). Planck should improve this 
constraint with its better measurement of the anisotropies on smaller scales; a detection 
limit of /nl ~ 5 is forecasted [163]. 

6. SECONDARY ANISOTROPIES 

Secondary anisotropies are generated after recombination. They can arise from grav- 
itational effects (the integrated Sachs-Wolfe effect or gravitational lensing) and from 
scattering events during or after the epoch of reionization. We have already met some 
examples, such as the large-angle polarization generated from scattering at reioniza- 
tion (Sec. 3.3). In this section we briefly review the other main sources of secondary 
anisotropy and their cosmological potential. For a recent, thorough review of the sub- 
ject, see [13]. Fully exploiting the primary anisotropies on multipoles / ~ 2000, for 
example to tighten constraints on Us, requires careful modelling of secondary effects. 
Such modelling is also important since the secondary signals themselves contain valu- 
able cosmological information about structure formation at late times and reionization. 
A number of ongoing ground-based surveys with arcmin resolution are seeking to ex- 
ploit the secondary anisotropies in this way [165, 166, 167, 168]. 



6.1. Gravitational secondary effects 

6.1.1. Integrated Sachs- Wolfe effect 

The integrated Sachs-Wolfe (ISW) effect is described by the last term on the right 
of Eq. (14). It is an additional source of anisotropy due to the temporal variation of 
the gravitational potentials along the line of sight. If a potential well were to deepen 
as a photon crossed it, the photon would receive a net redshift, and there would be 
a decrement to the CMB temperature along the line of sight. Conversely, a decaying 
potential well gives a temperature increment. At late times, the potential evolves only 
when dark energy (or curvature) start to dominate the background dynamics - a linear- 
theory effect - or as non-linear structures form. 

The linear effect from evolving potentials contributes to the temperature anisotropies 
only on large scales (see Fig. 11). This is because there is little power in the gravitational 
potentials today on scales that entered the Hubble radius during radiation domination, 
and, furthermore, as an integrated effect, the contribution of small-scale fluctuations is 
suppressed. The physics of dark energy affects the decay of the potential through the 
background expansion rate (via its energy density and equation of state) and through 
its clustering properties (i.e. sound speed). Both impact on the amplitude and detailed 
scale-dependence of the ISW contribution [169]. The late-time ISW effect is the only 
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FIGURE 11. Predicted secondary anisotropies from gravitational effects (left) and scattering effects 
(right) from [21]. (a) The ISW (blue), Rees-Sciama (purple) and lensing (red; negative where dashed) 
effects compared to the primary anisotropies from recombination (green) and the sum of primary and 
gravitational secondaries (black), (b) The sa —Tie fractional suppression of the primary anisotropies 
(red; discussed in Sec. 2.7.2), linear Doppler effect (blue), thermal Sunyaev-Zel'dovich effect (purple; 
calculated in the Rayleigh-Jeans limit), the density-modulated Doppler effect with linear densities (thin 
cyan) and the total including non-linear modulation e.g. from clusters of galaxies (thick cyan), and 
the Doppler effect modulated by the patchiness of reionization ("i-mod"; purple) compared with the 
primary anisotropies (green). The predictions for the SZ effect, non-linear Doppler-modula tion and patchy 
reionization are based on simplified halo models; see [21] for details. The quantity AT" = \Jl{l + l)C//2;r. 



way to probe late-time structure growth with linear CIVIB anisotropies alone, but this 
is hampered by cosmic variance on large angular scales. A more promising route is 
to cross-correlate the CIVIB with a tracer of large-scale structure [170]. A positive 
correlation is expected from decaying potentials, and was first detected by correlating 
the first-year WIVIAP data with the X-ray background and the distribution of radio 
galaxies [171]. Since the initial detection, the correlation has been detected with a 
number of large-scale structure tracers. The most complete analyses to date [172, 173] 
find ~ 4.5c7 detections. The power of the ISW effect in constraining dark energy models 
is rather limited by chance correlations between the primary CIVIB and large-scale 
structure; the total signal-to-noise can never exceed ~ 7, even for a tracer perfectly 
correlated with the ISW signal. Current detections provide independent evidence for 
dark energy, but no evidence for departures from ACDIVI (i.e. dynamical dark energy). 

The non-linear late-time effect is sometimes called the Rees-Sciama effect [14]. It 
arises from the non-linear growth of structures and, more generally, from the bulk motion 
of clustered matter. The power spectrum of the Rees-Sciama effect is broad but is a sub- 
dominant effect on all scales [174] ; see Fig. 1 1 . 
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6.1.2. Weak gravitational lensing 



Weak gravitational lensing of CMB photons by the large-scale matter distribution at 
z < 10 has several important effects on the CMB; see [17] for a recent review. The deflec- 
tion arises from the cumulative effect of large-scale structure along the line of sight. The 
total deflections are predicted to have an r.m.s. ~ 3 arcmin, and be coherent over scales 

1° (the angle subtended by the typical coherence size of the gravitational potential, 
~ lO^Mpc, at a comoving radial distance where lensing is most efficient). Lens- 
ing conserves brightness, simply re-mapping points according to the deflection field. 
For the temperature anisotropies, this re-mapping leads to a smoothing of the acoustic 
peaks resulting in fractional changes in Cj of around 10% at the troughs by / ~ 2000. 
Moreover, it transfers large-scale power to small scales so that it dominates the primary 
anisotropies for / > 4000; see Fig. 11. The effects in polarization are similar, but with 
the important addition that 5-mode polarization is generated from E [175]. This arises 
since the re-mapping by lensing does not preserve the geometrical relation between po- 
larization direction and the angular variation of the amplitude which defines E-mode 
polarization. The resultant B modes have an almost white spectrum for / <^ 1000 with 
Cf r^lx 10~^juK^; see Fig. 4. Characterising the 5-mode polarization from lensing is 
a secondary science goal of upcoming polarization experiments that have gravitational 
waves as their primary target [176]. 

As the angular resolution and sensitivities of CMB experiments have improved, it has 
become important to include the effect of lensing when deriving parameter constraints 
from CMB temperature data. For example, in the recent analysis of ACBAR data [30], 
a significantly better fit of ACDM models to the data is found when lensing is included. 
Including the amplitude of the lensing correction as a parameter in the analysis (with a 
predicted value q = 1), gives q = 1.23+|];^^ (95% C.L.). 

An additional feature of lensing is that it introduces non-Gaussian signatures in the 
CMB. Roughly, sub-degree scale CMB features behind a given lens are sheared coher- 
ently producing what locally looks like a statistically-anisotropic pattern of Gaussian 
CMB fluctuations. Averaged over lenses, statistical-isotropy is restored but the lensed 
CMB becomes non-Gaussian. This non-Gaussianity, manifesting itself as a non-zero 
trispectrum [32], may be used to reconstruct a noisy estimate of the deflection field 
(e.g. [177]). The noise arises from both instrument noise and chance alignments in the 
primary (random) CMB field that mimic the effect of a lens. The more small-scale fea- 
tures in the field being lensed, the better the reconstruction. Current data is too noisy to 
allow a direct detection of lensing via the reconstruction route. The signal-to-noise can 
be improved by cross -correlating the reconstruction with a tracer of large-scale structure. 
[For a high density of sources, the lower (Poisson) noise compensates for the reduced 
level of correlation.] The lensing effect in the CMB temperature has now been detected 
at the 3c7 level in cross-correlation with galaxy surveys [178, 179]. 

Moving beyond these first detections, the interest in CMB lensing will turn to exploit- 
ing it as a probe of the growth of structure; see [176] for a recent overview. In particular, 
CMB lensing has the potential to provide constraints on dark energy, modifications of 
gravity, and massive neutrinos that are inaccessible with the primary CMB anisotropies 
alone. CMB nicely compliments cosmic shear surveys [180] as it probes structure at 
higher redshift and with quite different systematic effects. 
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Note finally that since the effect of gravity on photons is achromatic, the frequency 
dependence of all gravitational secondary anisotropies follows that of the primary 
anisotropies. 



6.2. Scattering secondary effects 

6.2.1. Thermal Sunyaev-Zel 'dovich effect 

The thermal Sunyaev-Zel'dovich (SZ) effect arises from Compton up-scattering of 
CMB photons by electrons in hot, ionized gas such as in galaxy clusters [18]; see 
also [181] for a review. The CMB temperature at the redshift of the cluster is much 
less than the gas temperature, ^b^" ~ fewkeV, so photons gain energy, on average, by 
scattering. The net fractional change in the effective temperature of the CMB is 



0(n;v)=3;(n) 



hv \ / hv \ ^' 
coth — — -4 



TqmB J V TqmB J 



(107) 



where the Compton y parameter is determined by the integral of the gas pressure along 
the line of sight: y{n) = {oi/nig^c^) J p^dl where mg is the electron mass. Unlike the 
gravitational secondaries, the thermal SZ effect has a characteristic frequency depen- 
dence - negative below the null at 217 GHz and positive above. For a typical galaxy clus- 
ter, ~ 10""^ around 100 GHz giving a decrement of O(IOO) jUK. In the Rayleigh- Jeans 
limit, hv <^ ^bTcmb. the decrement is independent of frequency: 0(n; v) ~ —2y{n). 

The SZ effect is now routinely detected in targeted observations towards known 
clusters. A significant feature of the SZ effect is that the observed temperature change 
is independent of the redshift of the cluster making SZ also a promising route to search 
for high-redshift clusters. Indeed, a number of blind surveys are now underway, e.g. 
AMI [165], SZA [166], SPT [167] and ACT [168], to do just this. These ground-based 
instruments, covering a broad spread of frequencies, all have the high resolution (~ 
1 arcmin) required to resolve clusters, and, together, expect to find thousands of clusters. 
The Planck satellite, despite its poorer resolution, will produce a full-sky catalogue with 
0(1000) SZ clusters. Recently, the detection of three previously unknown clusters was 
reported by the SPT team [182]; these are the first clusters to be found in a blind SZ 
survey. 

Once detected, clusters will be foUowed-up with redshift measurements to determine 
cluster counts as a function of redshift, d^'N /dzd£l. This is related to the underlying 
cosmology through the expansion rate H{z), the cluster mass function, d«(M,z)/dM 
giving the number density of clusters per mass M as a function of redshift, and the 
survey selection function. The mass function depends sensitively on the growth of 
structure making cluster counts a promising route to improving our understanding of 
dark energy [183]. However, to achieve this goal, the survey selection function must 
be very well understood and this depends on the relation between the SZ observables 
and the cluster mass. Fortunately, the total SZ flux is expected to be a good proxy 
for the cluster mass, with the two being related by scaling relations. We can see that 
this is reasonable by noting that the SZ flux is proportional to the volume integral of 
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the electron pressure, by Eq. (107), and hence electron thermal energy. For a relaxed 
cluster, the virial theorem relates the thermal energy, hence SZ flux, to the virial mass. 
For a flux-limited SZ survey, the selection function should be rather uniform in mass. 
Better understanding the relation between mass and SZ observables, both with further 
simulation work and X-ray and optical follow-up, is an important aspect of SZ cluster 
research. 

The thermal-SZ effect is an important foreground for studies of the primary CMB 
temperature anisotropics on small scales. The power spectrum is shown in Fig. 11 in 
the Rayleigh- Jeans limit. On large scales it is essentially white noise, but turns over 
on scales where clusters are resolved. Since the frequency spectrum is known, and 
differs from the primary CMB anisotropics, the thermal SZ effect can be removed with 
multi-frequency observations or by observing close to the null at 217 GHz. Moreover, 
the signal is very non-Gaussian. However, for a survey like Planck, masking clusters 
detected at high significance in the CMB maps themselves is not sufficient since the 
residual contribution from high-redshift, low-mass clusters below the detection threshold 
is still sizeable [184]. The power spectrum of residual thermal SZ emission may already 
have been seen by CBI [185] and BIMA [186] as an excess of observed small-scale 
power over that expected from the primary anisotropies. However, to reconcile the 
observed excess with SZ requires around twice the number density of clusters than 
predicted for standard cosmological parameters. Interestingly, recent data at the same 
frequency (30 GHz) from the Sunyaev-Zel'dovich Array [187] find no evidence for 
excess power after correction for the contribution of residual point sources. The origin 
of the discrepancy between these results is still unclear. 



Thomson scattering of CMB photons off free electrons in bulk flows after reionization 
generate frequency-independent temperature anisotropies 



via the Doppler effect. For an optical depth, Tre ~ 0.1, and r.m.s. peculiar velocities 
~ 10~^, the linear Doppler effect would swap the primary anisotropies were it not for 
a rather precise cancellation from integrating along the line of sight. The cancellation 
follows since n ■ Vb is the radial gradient of a potential in linear theory. The small net 
effect is shown in Fig. 1 1 and is at most two orders of magnitude below the power 
spectrum of the primary anisotropies. 

Any spatial modulation of the electron density will tend to avoid the line-of-sight 
cancellation. This can arise from patchiness of the ionization fraction during reioniza- 
tion [188], or from the dumpiness of the (almost) fully-ionized intergalactic medium 
after reionization is complete [18, 20, 189]. Both effects produce anisotropies on scales 
of a few arcmin and below with the same thermal spectrum as the primary anisotropies. 
Computations of the power spectra of the modulated Doppler signals are shown in 
Fig. 11. The details of the patchy signal depend on the morphology of the ionized frac- 
tion (which affects where the power peaks) and the duration of reionization. Since the 
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6.2.2. Doppler effects 




(108) 



physics of reionization is still rather uncertain, and calculating the anisotropics from 
patchy reionization properly requires large numerical simulations incorporating radia- 
tive transfer for the ionizing radiation, there remain significant modelling uncertainties 
in the predicted power spectrum.^'' However, the patchy signal is expected to be sub- 
dominant to the density-modulated Doppler signal. There are contributions to the lat- 
ter from the linear density field (often called the Ostriker-Vishniac effect [20, 189]), 
which is easily calculated analytically, and from non-linear structures such as clusters of 
galaxies and filaments at low redshift (the kinetic Sunyaev-Zel'dovich effect [18]). The 
density-modulated Doppler signal is dominated by the non-linear effect on small scales 
(/ > 1000). 

Observations of the small-scale temperature anisotropics can provide complementary 
information on reionization to that from large-angle polarization. In particular, the latter 
is mostly sensitive to the integrated optical depth Tre, with only weak sensitivity to the 
duration of reionization [60], and none to the morphology. By having one observing 
channel close to the null of the thermal Sunyaev-Zel'dovich effect, experiments like the 
South-Pole Telescope and the Atacama Cosmology Telescope should detect the power 
spectrum of the Doppler signals at very high significance [190]. 

Small-angle polarization generated at reionization should be around two orders of 
magnitude below the temperature signal generated there and will likely remain unob- 
servable [191]. Note, however, that it is a further source of 5-mode polarization. 
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